naa Seller ele etch a pra ene rings ; . : Sdemeseeernnone : pheet ae 


For Reference 


= 
cS) 
) 
ee 
z 
é 
Be 
Z 
fal 
e- 
< 
a 
ty 
69 
g 
SS 
S) 
Z 


Gx wenis 
UNIOKASTTATIS 
LARERTAEDSIS 


SSS Sess 
The University of Alberta 
Edmonton, Alberta 
Printing Department 


ey Pay Der ay ' A i Ate os 
' - nS WAY 
Pee Ny { Ns Ci), ne G Re 
, ‘, th ‘7 
” ii V} 
ah j 
f 5 
iV ‘ 
oe a 
yea? * : ; = 
rh fi 
} 
i 
iii. | ry hoy) | on Ue 7 7 ey 
i he be, a WO ua * wat 
‘A i \ \ ae : Ee 
i : 7 ! 
is 
{ Wiz 
aay 
x. 
. 
” E 
ar | 
A . 4 
i) om. 
CP 
, ari 
ant ; 
7h j 
x wy \ 
t ar at tp ry 
Taps . 
Ny La \ - 
i : ir ; 
; : ; ys i AL 
: : : 
‘ - 4, Me r ‘ a » 
ri ange 
7 / oy ¥ 
" ry toe 7 “ 
_ 
i et 7 i 
i 
7 


THE UNIVERSITY OF ALBERTA 


PHYSICAL TAPERING OF LARGE ANTENNA ARRAYS 


by 


(C) pestows ANTHONY WYNNE 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIRMENTS FOR THE DEGREE 
OF MASTER OF SCIENCE 
IN 


ELECTRICAL ENGINEERING 


DEPARTMENT OF ELECTRICAL ENGINEERING 


EDMONTON, ALBERTA 


Peleg eS aes 


DEDICATION 


To my Parents 


iy 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Wynne1973 


ABSTRACT 


The performance of large fan beam arrays (160) by 
4) for T-type Radio Telescopes is studied using a summation 
technique that takes advantage of the fast Fourier transform 
algorithm. The effects, on sidelobe levels and beamwidths, 
of physical and resistive tapering of the excitation 
distributions are examined. The effects of random errors in 
the excitation of the elements are also considered. It is 
shown that the number of elements required for this type of 
array system may be substantially reduced, without degrading 
the performance of the system, by using a pow inah en of 


both physical and resistive tapering. 
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CHAPTER 1 


INTRODUCTION 


1.1 General 

In recent years there has been considerable 
interest in low frequency radio astronomy. Work has been 
done at 10 MHz and 22 MHz at Penticton in western Canada 
[1,2]. A broadband system from 10 MHz to 25 MHz is in 
operation in central Russia [3,4]. ‘The University of 
Maryland's Clark Lake Radio Observatory has a system 
presently under construction which will operate between 10 
MHZ and 110 MHz [5]s however, below 20 MHz its use will be 
limited as the collecting area of the antenna will be low. 


In the southern hemisphere work has been done at 3.5 MHz and 


ae hhz. [ 6,7). 


All of the above radio telescopes use a cross or 
tee type antenna system first described by Mills and Little 
{8}. This type of telescope produces two orthogonal fan 
beams and by the multiplication of these beams a pencil beam 


is obtained which gives high resolution at low frequencies. 


At present there is a need for a new radio 
telescope that has greater resolution and sensitivity at low 
frequencies than any of the instruments presently available. 
A proposal for such an instrument, to be located in Alberta, 
has been made and site testing is being carried out [9]. 


The telescope would operate near 12 MHz and would consist of 
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an array of horizontal antennas arranged in a tee with a 
north-south arm 2.5km by 100m and an east-west arm 5km by 
50m. Such an array would have a collecting area of the 
order of 100,000m2 and a resolving power of about 30 arcmins 


at the zenith. 


The distribution of radio emission in the sky is 
very complex. Radiation is received from sources of all 
angular sizes in addition to the general background 
emission. For this reason a pencil beam antenna system with 
as low a sidelobe level as possible is most desireable. 
Sidelobe level reduction is achieved through what is 
referred to as tapering [10]. ‘Tapering consists of 
gradually reducing the level of the excitation of the array 


elements to some small value at the ends of the array. 


There are many forms that the excitation 
distribution may take. A uniform excitation distribution 
produces a narrow beam but the sidelobe levels are high. 
For large arrays the level of the first sidelobe is 13.6 dB 
below the peak of the main beam. A beam with no sidelobes 
could be produced by using a binomial excitation 
distribution but this increases the beamwidth and reduces 
the gain of the array by a large amount. Other excitation 
distributions such as triangular, cosine or cosine squared 


produce beams with intermediate characteristics. 
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In the following chapters, array systems which are 
160 wavelengths long, less than 4 wavelengths wide, and that 
have an excitation distribution approximating a Gaussian 
will be examined. Such arrays produce fan beams in a 
direction transverse to the principle axis of the array. A 
true Gaussian excitation distribution, infinite in extent, 
would produce a beam with no sidelobes. In practice, the 
excitation distribution of a finite array can be represented 
by a truncated Gaussian. Choosing a taper such that the 
excitation of the elements at the ends of the array is 1/16 
that of the centre elements [11] represents a good 
compromise between the requirements for gain, beamwidth and 
Sidelobe levels of the power pattern. Some typical 
excitation distributions are illustrated in Figure 1.1 for a 
linear array 160 wavelengths long. Note that the excitation 
is assumed to be symmetric about the centre of the array. 
In each case the excitaticns of the centre elements have 


been set to unity. 


The most straightforward method of obtaining a 
specified excitation distribution is to place attenuators in 
the lead to each antenna element. This approach is known as 
resistive tapering. Another technique is to completely 
eliminate some of the elements of the array or to replace a 
group of elements, in a Section of the array where the 
excitation is low, by a smaller group of elements with 


modified excitations. This is known as physical tapering. 
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The excitations of the elements in physically tapered arrays 
may be found by setting the total contribution of the 
elements in sections of the array egual to the total 
contribution of the elements in the same sections of a 
resistively tapered array. It is also possible to use a 
combination of both physical and resistive tapering in an 


array to achieve a specified excitation distribution. 


In the case of large array systems, physical 
tapering offers many practical advantages over resistive 
tapering. The physically tapered arrays discussed in 
Chapter 4 have 37% fewer elements than the corresponding 
resistively tapered arrays. This reduction represents a 
Significant saving in the cost of constructing such an array 
system. If the system requires signal processing at each 
element there will be a further saving in the processing 
equipment needed. The time involved in setting up and 
tuning the array will also be reduced. There will be fewer 
points at which errors can occur and if everything else is 
unchanged the effect of the errors should be less on a 


physically tapered array. 


1.4 Purpose of the Thesis 

The purpose of this thesis is to compare the 
performance of resistively and physically tapered arrays in 
which the widths of the arrays are small compared to their 
lengths. A technigue will be developed in Chapter 2 for 


finding the power pattern of one and twce-dimensional arrays. 
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It will be shown that if all the elements in the array are 
identical then the power pattern can be split into the power 
pattern of a reference antenna times an array factor which 
takes into account the location and excitation of every 
elememt in the array. A comparison between different array 
systems can then be made on the basis of their array 
factors. In Chapter 3 the effects of coarse resistive 
tapering on the gain, beamwidth, average sideiocbe level and 
peak sidelobe levels of one-dimensional arrays will be 
examined. These results will be extended to two-dimensional 
arrays, in Chapter 4, and a comparison between resistive and 
physical tapering will be made. The difference between 
resistive and physical tapering will be shown to be greatest 
at points far from the main beam. The effects of random 
errors in the excitation of the elements of both one and 
two-dimensional arrays will be studied in Chapter 5. It 
will be demonstrated that the effects of errors are greatest 
on the carefully adjusted resistively tapered arrays and not 


so great on the physically tapered arrays. 
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CHAPTER 2 


ANALYSIS OF THE POWER PATTERN 


The power pattern of an array of antennas may be 
determined by summing the contributions from each element at 
each angle for which the total field is to be calculated. 
This method is not practical for large arrays because of the 
time involved in doing the calculaticn and the problems 
encountered in retaining accuracy. Considerable 
Simplification may be achieved by use of the concept of the 
array factor which is the power pattern of an array of 
isotropic point sources with the same physical arrangement 
as the actual array. The following analysis of the power 
pattern and the array factor for rectangular arrays is 
directed towards finding the result in a suitable form for 
calculation using computing facilities such as those 


provided by the IBM model 360/67 computer. 


Throughout the analysis the effect of the mutual 
impedances between the elements will be neglected and it 
Will be assumed that each of the elements in the array can 


be excited as specified. 


The analysis of the arrays will be performed using 
Fourier transforms and so the behaviour of the arrays will 
be examined in the transform space (K-sSpace). The primary 


advantage of working in K-space is that the fast Fourier 
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transform algorithm [12], which is both fast and accurate, 
may be used to compute the array factor. Another advantage 
is seen when dealing with phased arrays. An unphased arrary 
is one in which the phases of the excitations of every 
element are identical. Beam steering in K—space is 
represented by simple linear shifts as will be shown in 
Section 2.8. Although the array factor can be transformed 
back to real-space, the process would involve a lot of 
computing time and for the purpose cf comparisons between 


different array systems it is not really necessary. 


The analysis of arrays with physical tapering will 
be performed by breaking the arrays up into sub-arrays where 
each sub-array is of constant width. The total field 
strength pattern may then be found by summing the field 
strength patterns of the sub-arrays. For array systems in 
which the widths of the arrays are very much less than their 
lengths this technique achieves an additional saving in 
computational time over a direct application of the fast 


Fourier transform. 


The calculation of the array factors will involve 
the generation of large tables of data and a repetitive 
sequence of operations at about 106 points in K-space. The 
computational time can ke reduced considerably by efficient 
coding of the repetitive sequence in assembler 


language { 13]. 
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2.2 Radiation 


Ss 


Radiation in unbounded space from a local source 


may be described in terms of a magnetic vector potential, A, 


which satisfies the Helmholtz equation 2-1 [14]. 


—_ 


-J 


V2A + k2Aa (2-1) 


J is the current-density (Amps/m?2) and k is the wave number 


of the medium and is given by 
(2-2) 


The resulting magnetic and electric fields (E and H) are 


functions of A as indicated in equations 2-3a and 2-3b. 


E = -jwuA + _1. V(V-A) (2-3a) 
J we 
H = VXA (2-3b) 
A(z) = 1 f f f d(E*) Bxp{-jk{T-E' [pax'ay'az' (2-4) 
Ut co -c -0 itor s lana aA 


Equation 2-4 is a general solution of equation 2-1 for A in 


an unbounded region. The primed variables represent the 


source coordinates and the unprimed variables the field 


coordinates. 
coordinates. 


localized are 


The integration is carried out over the source 


When the current distribution in some 


a is known, 


equations 2- 


3a, 2-3b and 2-4 


completely specify the resulting electromagnetic radiation. 
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2.3 Radiation from a Y-directed A/2 Dipole in the X-Y¥ Plane 
Consider a Y-directed \/2 dipole located at X',yY' 

in the X-Y plane as indicated in Figure 2.1. Assuming that 

the dipole is oriented with its axis parallel to the Y-axis, 

the current distribution and hence the magnetic vector 

potential of equation 2-4 will have only a Y-component. 

When only the nature of the field at large distances from 


the antenna is required, |{f-r'| may be replaced by 


(Appendix A) 


> —> 


{r-rc'| = r{1-Sin6d (X'Cos@+Y'SinQ@) } (2-5) 
rc 


Thus equation 2-4 becomes 


J J J dye") Exp (-jkr+jkSine (x! Cosé+¥' Sing) } dx" dy'dz 


Ay(r) = 1 
r-Sind (X 'Cos@t+Y'Sing) 


1 
(2-6) 


In the denominator of equation 2-6 the trigonometric terms 
are small compared to r and may be neglected. In the 
exponent they represent phase shifts and since X' and Y' may 


be comparable to \ they must be retained. 


Let Kx, = kSin8Cos¢ (2-—Ta) 
hg = kSindSing (2-7b) 
A(r) = Exp {-jkr} (2-8) 

Gnr 


Then equation 2-6 becomes 
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Figure 2.1. | Coomdinatej System: 
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Figure 2.2. Location of the Elements. 
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Spherical Coordinate 6 


The relationship between the parameters of K-space and real 
Space are shown in Figure 2.3. kK, and Ky extend to infinity 
in all directions but, as will be shown in Section 2.6, only 


a limited area of K-space represents real space. 


Let the integral part of eguation 2-9 be denoted 


by 


A(KxeK 


pee | i dy eine ieee de ay faz" (2-10) 


Then Ay (Lr) becomes 
Ay (CT) = A(r)A (Ky,/Ky) (2-11) 


In equation 2-10 Jy (r‘) is a current density in Amps/m2. 
Consider the dipole, of length L, located at (X!$,Y¥{) in the 
X-Y plane as indicated in Figure 2.2. The subscripts m and 
n are introduced at this point for convenience. In Section 
Zoo (48,15) will refer toe the position of the element \(#,n) 
in an array of elements. If I(X{$,Y§) represents the total 
current flowing along the dipole then Jy (r") may be written 


as 
Jy(E") = T(Xhe YA) 6 (X*- Xp) 9 (Z") (2-12) 


where the delta functions have the dimensions of inverse 
length. It may be assumed that the current at the ends of 
the dipole is zero and varies sinusoidally over its length. 


Hence, Jy (t') becomes 


Sy (E%) = IoCynd (XI MEU (NATE) (297 Sin (kr) (2-13) 
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where I, is the peak current on a hypothetical antenna and 


where 
Cnn = the current on a dipole located at (Xi,Y3) 
To out ad yeera(2s4s) 
U(Y'-¥Y#) = u(¥'-Y"-L/2) - u(¥'-Y4+L/2) (2-15) 
Y= L/2-|Y'-y"] (2-16) 


Cnn has been introduced at this stage so that in subsequent 
sections arrays of dipoles with different currents may be 
discussed directly in terms of the results presented here. 
For this current element located at Xj}, and Y{ the magnetic 


vector potential becomes 


Aye (= EGY Ainy (Rye Ky) (2-17) 


where Ann (Ky Ky) is obtained by substituting equation 2-13 


into equation 2-10. Thus 


CS ate Oe Ke ee mno (X'-X8) U(¥'-¥4) 6 (Z") Sin(kYy 


Exp{jK,X'+jK,¥"} JdXtayeaz' (2-18) 


The indicated integration is carried out in Appendix B where 


it is shown that 
Bian (Ky eKy) = Cop EXP (5KyX gti Ky Va} Ayr (Kye Ky) (2-19) 
where 


an CEA ek - Cos {kL/2} ] (2-20) 
Sin@@ Sine? 
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Ayr (Kys Ky) is known as the element pattern. It may be 
interpreted as the magnetic vector potential that would 
result if the element at (At,-Y4) were located at the centre 
of the coordinate system with its peak current adjusted to a 
hypothetical reference level, I,. The exponential factor in 
equation 2-19 represents the phase shift associated with 
displacing the element from the origin to (X!,¥'). The 


constant C,, accounts for the relative magnitude and phase 


n 
of the current in the element (m,n) with respect to the 


reference level, Tye 


It was assumed in the earlier sections that the 
dipole was radiating in free space. In practice the dipoles 
will be placed over a reflecting screen which will modify 
the fields and power pattern of the antenna. These effects 
will be considered separately as it will be shown that they 


are common to all the different array systems to be examined 


The effects of the screen may be calculated by 
considering the antenna as composed of two antennas one 
above the screen and the other, an image antenna with its 
current reversed, below the screen. The reversal of the 
currents is necessary to produce the zero tangential E-field 


which must occur on the screen. 
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For the antenna above the plane of the screen the 
phase of the magnetic vector potential must be advanced by 


an amount YW where 
vw = khCosé@ (2-21) 


and h is the height of the antenna above the screen. For 
the array below the screen the phase must be retarded by the 
same amount. By combining the two antennas eguation 2-17 


becomes 
Aymn (T) = A (IL) Amn (Ky Ky) [Exp {jk hCos 6} -Exp {-jkhCos6 ] (2-22) 


This additional factor will be called the screen factor 


SF(K,) which may be written as 


SF(K,) = 2jSin (jhK,) where K, = <xCos0 (2-23) 
Thus, 
Agon(o) = Uden ABE (K, Ky) SF (Kz) (2-24) 


SF(K,) depends only on the distance from the reflecting 


screen to the antenma and not on the shape of the antenna. 


2.5 Radiation from an Array of \/2 Dipoles in the X-¥ Plane 
Assume that the array is rectangular and that each 
position in the array contains a dipole. This assumption 
produces no loss in generality, since, for a non-rectangular 
shaped array it is possible to fill out the array toa 
rectangular shape by including a number of additional | 
elements for which the level of the excitation is set to 


zero. The resulting fields and power pattern will be 
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unchanged. Assume also that the numbers of dipoles in the X 
and Y directions are even. For the case in which the number 
of dipoles in either direction is odd only a slight 
modification to the method of analysis described in this 


section and the following sections is required. 


Let the array of dipoles be laid out as shown in 
Figure 2.4. The spacings between the dipoles in the X and Y 
directions are D, and Dy respectively and the numbers of 
elements in the X and Y directions are N, and Ny 
respectively. Let (m,n) be the system used for identifying 


the element located in column m and row n. (XP, YE) isthe 


location of element (m,n). 


Ay () for the complete array is found by summing 


By mn (©) over all m and n. Thus from equation 2-24 
Nx-1 Ny-1 
9 ea) ; A (2) Ann (Ke Ry) SF (K2) (2-25) 
m=O n=0 
Nx-1 Ny-1 
= A(r)SF(K,) J y | Ang(@e Sy) (2-26) 
m=O n=0 
From eguation 2-19 
Nx-1 Ny-1 
Ay(r) = A(x) SF(Kz) Ayr (Kx Ky) ) ) CmnExXp{jKxXhtjKyY h} 
m=0 n=0 2720) 


By defining an array facfor AF(K,,Ky) as 


Nx-1 N -1 ; 
as ) ConExP (gh xXhtahyi ht (2-28) 


AF (K,,K 
m=O n=0 
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Figure 2.4. Arrangement of the Elements in an Array. 
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the magnetic vector potential for the complete array becomes 
Ay (©) = A(L)Ayr (Ky Ky) SF (Kz) AF (Ky, Ky) (2-29) 


Equation 2-29 has four distinct terms. A(r) represents the 
dependance on the distance from the antenna to a point in 
the far field and will nct be taken into account in the 
comparison of the different array systems. The element 
pattern Ayr (Ky eKy) - takes into account the shape of the 
elements used in the array and the basic form of the current 
distributions in the elements. The screen factor, SF(K,), 
takes account of the effects of a reflecting screen placed 


below the antenna. The array factor, AF(K,,K accounts 


ae 
for the physical layout and arrangement of the elements in 
the array. It also takes into account the differences in 


the magnitude and phase of the excitations between the 


elements. 


The E and H vectors now can be found from 


equations 2-3a and 2-3b. 


Ey = ~JWOUA, (2-30a) 
ee) sig Mea 2—=30D 
Paced 4 ) 
ae 23 
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Ag = Ay(r) Cos@Sing (2-3 1a) 
A = Ay (L) Cos (2-31b) 
Eg = ~jwuay (r) Cosésing (2-32a) 
Ey = ~jWuUA, (Lr) Cos¢ (2-32b) 
f,<> jka, (r) Cos@ (2-32c) 
Hy = ~jkAY (x) Cos@Sing (2-32d) 


Finally, the power pattern is given by 


uf * 
Br = (BQH Ey Ha) 


= [Ay(r)|24Uk(Cos29 Sin2@ + Cos2p) (2-33) 


The array factor is the only term which will vary 
when the arrangement of the array or the excitations of the 
elements are changed. Hence, comparison between different 
array systems will be made on the basis of their array 


factors. 


It will not be necessary to refer to the element 
pattern or the screen factor except when the array factor 
has large sidelobes near the zeros of the element pattern or 
the screen factor. Under such circumstances the element 
pattern or the screen factor may reduce these sidelobes to 
an acceptable level. Equation 2-33, for the power, contains 


an additional factor arising from the conversion of the 
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magnetic vector potential from rectangular to polar 

coordinates. This factor also contains zeros and if the 
array factor has large sidelobes at the location of these 
zeros this term must be taken into account in determining 


the actual sidelobe levels of the array. 


2.6 The Array Factor 
The location of the elements shown in Figure 2.4 


is given by 


X' = (m-Nx-1) Dx (2- 34a) 
2 

Y§ = (n-Ny-1) Dy (2-34b) 
2 


By substituting these values for X/ and Y! in equation 2-28 


the array factor becomes 


Nx7-1 N -1 
AF (K,,K,) = BKK) ) ; Con EXP{jK,D,m+jK Dyn} (2-35) 
m=O n=0 
Where 
B(K,/K,) = Exp {-jK,D, (N,-1)/2 -JjK,D, (Ny-1) 72} (2-36) 
27 
Py = KyD (2-37b) 
“or 
Then, 
Nx-1 N = P 
AF (K,/K,) = B(K,,K,) : j CanEXP{j27P,m+j2mP pn} (2-3 8) 
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The advantage of expressing the array factor in 
this way is that the summation part of equation 2-38 is very 
Similar in form to the two-dimensional finite Fourier 
transform [15] which may be written as 

Ni — 1 No -1 
X(Liel2) =  ) =) A(Ky,K2) EXp{52™KyL,/N, +527K,L2/N3} 
K,=0 K5=6 (2-39) 
Where N,, and N, are integer constants and where L,, Io, Ky, 
and K, are integer variables defined in the following 
ranges. 


@.<5 Ky Ss Ny-1 (2-40) 


oS Be S Ny-1 


X(L,,L,) and A(K,,K,) are complex variables. The summation 
defined by equation 2-39 may be computed using the 


techniques of the fast Fourier transform. 


Let | ee (2-41a) 
Ny = N, (2-41b) 

m= K, (2-41c) 

n= K, (2-4 4d) 

BY SOL 4/ny (2-41e) 

Boe LyN (2-41£) 

Can ‘ A(K,,K,) (2-419) 
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By using the substitutions defined by equations 2-41a 
through 2-41g it is possible to use the finite Fourier 
transform algorithm to do the summation part of equation 


2-38. 


It can readily be shown from equations 2-7a, 2-7b, 
22=37a and 2-37b that Kx, Ky, Px and Py are continuous 


variables which are defined in the following ranges 


-2n < Kx < 2n (2-42) 
r d 
s2m 5 Kyess 2m 
r r 
r r 
ee hy = Py 
r r 


It should be noted that although K,, Ky, Px and Py 
are continuous variables equations 2-41le and 2-41f define 
values of Px and Py which can be expressed only as rational 
numbers. For this reason it is apparent that the array 
factor can only be found at discrete points in K-space by 
using the finite Fourier transform. In addition to these 
restrictions on Ky, Ky, Px and Py only certain values of 
these variables represent real points in the sky. These 


values may be found as follows 
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Sin@ = (K2+K2) V2 (2-43) 
Therefore 

o ee (heats) a. <ike (2-44) 
O < {(21Px)2 + (27Py)2} < k2 (2-45) 

ayes 7 D2 

x ¥ 
ré DY SAD SNt/e (2-46) 
then 0 < (Erte s) < 174 (2-47) 


The relationship between 6,6 and ehh is shown in 
Figure 2.3. The sidelobes of the array factor which occur 
in the regions of K-space outside the region defined by 
equation 2-43 do not contribute any power to the system. If 
the array is phased by the introduction of a progressive 
phase shift in the excitation of the elements, these 
Sidelobes may be shifted into the region of K-space 
representing real points in the sky and will now contribute 


power to the system. 


By uSing the technigues of the fast Fourier 
transform the array factor can only be found at points in 
K-space where P, and Py are rational numbers as was 
indicated in the last section. The intervals in P.. and P 
between the points ded thee by equations 2-41e and 2-41f are 
1/N; and 1/N2 respectively. Increasing Ni and Nz reduces 
the size of these intervals and increases the resolution 


since the total range of P, and Py is limited by equation 
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2-47. WN, and N, will both be set to 2048 in the following 
work. This choice sets the increment in K-space to 


0.006136/) for both K, and K, when D,=Dy=\/2. Also, in this 


be 
case, it sets the overall size of the array to 1024) by 
1024 which is much larger than any of the arrays to be 
examined in later chapters. The size of these arrays will 
be 160 by 4A. The array factor of these arrays can be 
examined with a high degree of resolution in K-space and 
hence in real space by using the fast Fourier transform on 


an array 1024 by 1024 where the level of excitation 


outside an area 160A by 4A is set to zero. 


Figure 2.5 illustrates how an array 320 elements 
long (160A) by 8 elements wide (4A) might be numbered. The 
excitation of the elements outside the shaded area is set to 


zero. In other words, 


Cc = 0 for-m<1020 and alin (2-48) 
for m>1027 and all a 
for n< 864 and all on 


for*n>1183 and alin 


For a physically tapered array the excitation of some of the 
elements inside the shaded area would also be set to zero as 
suggested by Figure 2-6. If the excitation of the elements 
does not change across the width of the array the n- 

subscript in C,, may be dropped and the array factor may be 


written as 
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Figure 2.5. Arrangement of the Elements in Resistively 
Tapered Arrays. 
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Figure 2.6. Arrangement of the Elements in Physically 
Tapered Arrays. 
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| if mo ‘ . cme : 


1027. 1183 
RE (ee eKy) = BK, Ky) } : CyExp{j2n (Pym+Pyn) } (2-49) 
m=1020 n=864 


By substituting for B(K,,Ky) and separating the summations 


with respect to m and n equation 2—49 becomes 


1027 
EP{k. eh) = G(P,) Bxp{-j2nP, (2047/2)3 ) Exp{j27P,m} 
m= 1020 (2-50) 
Where 
1183 
G(Py) = Exp{-j2mP, (2047/2) } ) ChExp{j2mP yn} (2-51) 
n=864 


In the actual running of the fast Fourier transform 
subroutine n will be incrimented from 0 to 2047, of course, 
since C,=0 for n<864 and for n>1183 the value of equation 


2-49 will be unchanged. 


AF (Ky /Ky) = G (PY) CExp{-j/7P,}+Exp {-j5mP,} ..-+Exp{+j/ TP} J 
= G (P,) 2{Cos (7mP,) +Cos (STP,) +Cos(37Py) +Cos (Px) } 
= (F(R {76(P)) (2-52) 
where 
F(P.) = 2{Cos (77P,) tCos (57P,,) +Cos (37P,) +Cos (™P,)} (2-53) 


It is apparent from equation 2-52 that when the 
width of the array is constant over its whole length the 
array factor may be expressed as the product of two 
orthoganal array factors. To find the array factor of an 
array which is physically tapered the array must be broken 


up into sub arrays in which the widths of the sub-arrays are 
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constant over their lengths. In the case of the physically 
tapered array shown in Figure 2.6 there are four distinct 
sections so there will be four terms in the array factor. 


Let the functions F(P,) and G(Py) be defined as follows. 


Ge({Py) = G(Py) where C, = 0 when n< 973 (2-54a) 
or n>1074 
Ge(Py) = G(Py) where Cy, = 0 when n< 945 (2-54b) 
or n>1102 
or 973sn<s1074 
Gy (Py) = G(Py) where Cy, = 0 when naon2 (2-54c) 
or n>1135 
or 945<n<1102 
G2kPy)= 6 (Py) i where ()Cqi=*0! when n< 864 (2-574) 
or n> 1is3 
OE 9t2<n24135 
Fe (®x) = F (Px) (2-55 a) 
Fe (Px) = 2{Cos (51Px) tCos(37Py) +Cos (TP) } (2-55b) 
F, (Px) = 2{Cos (37P,) +Cos (Px) } (2-55c) 
Fo (Pa) = 2 {Cos (TP, )) (2-554) 


The array factor for the array shown in Figure 2.6 can now 


ke written as 


AF (KysKy) = Fe (Px) Ge (Py) tFe (Px) Ge (Py) + Fu (Px) Gy (Py) +F2 (Px) G2 (Py) 


(2-56) 
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2.8 Beam Steering 

Beam steering may be achieved by progressively 
delaying (i.e. Phase Shifting) the excitation of the 
elements towards the desired position of the beam. The 
effect is to shift the array factor in K-space so that the 
Main beam points in the new direction. To show this effect 
consider beam steering along the principle axis. In this 


case @ is 90 degrees and equation 2-38 may be re-written as 


follows. 


Ny-1 
AF(KY)°= Exp {-jKyDy (Ny-1) 72} " C, Exp {jk Dyn} (2-57) 
n= 


For the case in which the phase of the excitation is zero 


set 


C STA (2-58) 


Ny-1 
AF (K,) Exp (-JKyDy (Ny 1) /2} J A, Exp {jK,D n} (2-59) 


For the case in which there is a progresSive phase shift in 


the excitation it is possible to represent C, as 


oO 
i) 


A, Exp {-JK, Dyn} (2-60) 


= 
in 
1) 
ue 
) 
w 
" 


kSin®, 


Ny-1 
2 pe -#h 22 A_E Ks K D 
AF(K,) = Exp{-jkyDy(Ny-1)/ ~ naa OURS B67 Saalpsetate 


If the array factor for the array with the progressive phase 


shift is now evaluated at Ky tk, instead of at Kye equation 
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2-61 becomes 


Ny-1 
bf 
AF (K_+K = Exp{-j(K_+K Dy (N24 j 
(RO tk! pirat Fagus y Ny ye dag eos 
Equation 2-62 may be written as 
AF(K_ +K = EX jJK oD €N =4) 72} shape (he — Ti 4Z 
eS.) PEK, (dD, (N,- 1) /2} Exp {-JKD. (N-1) 72} 
N= | 
A_E jK_D 2-63 
pis aExp (4K, Dn} (2-63) 


Except for the first exponential factor equations 2-63 and 
2-59 are identical. In other words, when the array factor 
for the case in which a progressive phase shift is present 
is evaluated at Ky +Ky, and the array factor for the case in 


which no progressive phase shift is present is evaluated at 


K the two array factors are identical except for a shift 


y! 
in phase. By taking the magnitude of both equations 2-59 
and 2-63 it is apparent that the shape of the array factors 
for the phased and the unphased arrays are identical and 
that the effect of introducing a progressive phase shift has 
been to shift the array factor in K-sSpace by an amount Kyoe 


Hence 


[AF (K,+K,.) | = [AF(KY) | (2-64) 


y°" Phased array Unphased array 

It should be noted that although the peak levels 
of the array factor do not change when the array is scanned, 
the widths of the main beam and the sidelobes do change 


which may produce changes in the relative power contributed 
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by the main beam and the sidelobes. The changes in width 
are due to the way 9 and @ are related to K, and Ky. For 
values of 6 close to zero in Figure 2.3 the smaliest 

increment in K, or Ky represents 0.056 degrees in 6 while 


for @ close to 90 degrees the same increment in K.. or Ky 


represents 3.532 degrees in @. It shouid also be noted that 


the peak levels of the power pattern may change when the 


array is phased. If this occurs it is due to the fact that 


the combined effects of the element pattern, the screen 
factor and the factor introduced by the change in 


coordinates from rectangular to polar do not produce an 


isotropic pattern in the direction in which the main beam is 


shifted. 
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CHAPTER 3 


In this chapter the array factor of five one- 
dimensional arrays will be examined. In all cases the array 
will consist of half-wavelength dipoles and their 
arrangement will be as indicated in Figure 3.1. 
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Figure 3.1. Arrangement of the Elements in One-dimensional 
Arrays. 
The tapering will ke entirely resistive and the excitation 
distributions will vary from uniform through four 
progressively better approximations to a Gaussian. The 
first three approximations to the Gaussian excitation 
distribution represent a very coarse taper. This type of 
taper will necessarily exist along the transverse axis of 
the physically tapered two-dimensional arrays to be 
discussed in Chapter 4. Indeed the exact form of the 
excitation distribution, in this chapter, will be chosen 
such that the response of these one-dimensional arrays will 
correspond very closely in their principle plane to the 
array factor of one of the four arrays to be discussed in 


Chapter 4. Correspondence is achieved by setting the 
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excitation at each point along the axis of the one- 
dimensional array equal to the sum of the excitations across 
the width of the two-dimensional arrays at the same points 


along the principle axis of these arrays. 


One-dimensional analysis allows the principle 
features of the array factor to be examined without the 
necessity of expensive computing. The beam width in the 
principle plane, the effects of phasing on the beam width, 
the location of the first zeros and the maximum sidelobe 


level that may be expected can easily be found. 


The analysis was done using an IBM 360 computer. 
Programs were written to generate the excitation data for 
the different types of arrays to be examined. ‘the summation 
part of the calculation was performed using a fast Fourier 
transform subroutine called HARM which is available in the 
IBM Scientific Subroutine Package, [16]. Additional 
programs were written to analyse the results and present the 
output in a useful form. These programs and the problems 


encountered in computing are discussed in Appendix C. 


3.2 Excitation 

The excitation data for the uniform and the 
truncated Gaussian arrays are shown in Figure 3.2a. The 
excitation is symmetric so only half has been shown. The 
shape of the Gaussian distribution is based on a taper of 
—24.0 dB over a distance of 80 wavelengths. The choice of 


—24.0 dB for the taper is based on a compromise between the 
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@ Uniform Array 


@ Gaussian Taper 


Excitation Level 
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Distance from the Centre ()) 


Figure 3.2a. EXcitation,Data fer la Unitorm and 
a Truncated Gaussian Array. 


@ 2-Stage Taper 


@ Gaussian Taper 


Excitation Level 
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Distance from the Centre (A) 


Figuré 3:2b: BxGitation bata for an Array with 
2 Stages of Tapering. 
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Figure’ 3.20. "Excitation vata Lor any ALray with 
4 Stages of Tapering. 


@ 8-Stage Taper 


@ Gaussian Taper 
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PiGure 3.20. ExXxCLeacion Data for an Array with 
8 Stages of Tapering. 
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requirements for gain and the maximum sidelobe level that 
can be tolerated. Figures 3.2b, 3.2c and 3.2d show the 
progressively better approximations from the uniform to the 
Gaussian excitation distribution. They represent tapering 
in two, four and eight stages where the excitation step 


sizes are 0.5, 0.25 and 0.125 respectively. 


In all cases the excitation of the centre elements 
has been set to unity so that the gain of the different 
systems may be compared. All the excitations are real. No 
progressive phase shifts were introduced as it was shown in 
Section 2.8 that the effects of Scanning on the array factor 


can be examined by shifting the array factor in K-space. 


3.3 Main Beam Response 

The array factors for the five array systems 
described in the last section were computed. Each was 
normalized with respect to its own maximum level. The 
normalizing factors are related to the relative gain between 


the systems. 


The main beam responses for levels greater than 
—6.0 dB are shown in Figure 3.3. It is apparent that the 
principle effect of tapering has been to increase the 
beamwidth. The response of an array with a true Gaussian 
excitation distribution was calculated but the response in 
this part of the main beam is not distinguishable from the 


response of the array with the truncated Gaussian excitation 


distribution. 
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Figures 3.4a through 3.4e show the first two 
degrees in the principle plane of each array factor. It can 
be seen that the beamwidth to the first minimum increases 
and the peak sidelobe level decreases as the number of 


Stages of tapering is increased. 


The response of an array with a true Gaussian 
excitation distribution is also shown in Figure 3.4e for 
comparison with the truncated Gaussian array. For 9 less 
than 0.8 degrees the difference between the two array 
factors is small, however, for 6 greater than 0.8 degrees 
the response of the true Gaussian falls rapidly to zero 
while the response of the truncated Gaussian has sidelobes. 
These sidelobes represent the effect of truncating the true 
GausSian excitation distribution and could be considered as 
representing the array factor produced by an error function, 
as suggested by Figure 3.5. The array factor of an array 
with type 1 excitation, as shown in Figure 3.5, is identical 
to the combined array factors of arrays with type 2 and type 


3 excitation. 


In all cases, for the array factors shown in 
Figures 3.4b° through 3.4d, for the arrays with coarse 
tapering, the locations of maximum sidelobe levels close to 
the main beam are not very predictable. The first sidelobe 
is not the greatest and the relative magnitudes of the 
sidelobes are not determined by the amount of tapering. 


There is, however, a general trend in the average sidelobe 
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Figure 3.5. 
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Excitation Data for Gaussian Arrays. 
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level which decreases as the amount of tapering is 
increased. The principle features of the main beam area are 
Summarized in Table 3.1. It is also instructive to define a 
‘figure of merit' for each array as indicated in equation 
3-1. Although this term will be referred to as the gain it 
Should be noted that it is not the directive gain usually 


associated with an antenna system. 


Gain = 20log,, AF (3-1) 
AFnax (Uniform) 


Table 3.1 
Characteristics of the Main Bean 


{| Taper { Gain | 1st Min [| Beamwidth | Max Sidelobe | 
| | (dB) | (Degrees) | (Degrees) | (dB) { 
| { { { i | 
a a a | 
iY pestorn |.’ 0.00%. 0.34 | 0.320 | a SIA, { 
| i i i { | 
| 2-stage | -2.34 .| 0.50 | 0.360 | -19.1 | 
| | { | { i 
{ 4-stage |{ -4.24 | 0.43 | 0.420 { -27.6 | 
{ | | { { i 
{ 8-stage | -4.84 | On 73 | 0.442 { Fae ARS, | 
i { | { i { 
{| Gaussian | -5.72 | 0. 84 | 0.468 | = 3-3 20 | 
| { | l | | 


The effects of phasing on the beamwidth were 
examined for each system by shifting the main beams, shown 
in Figure 3.3, to different locations in K-space and finding 
the new values of 9 which correspond to the 3 dB points 


located in the new diagram. The effects of phasing on the 


2h Daa THy sf 


SIs souk Waod ae 


p£ sukteb. oF ovttouisas: orks" i 


ot teipe ue o 


Fi othe. aAt, 26 


oe feet rr meee aemenmns slants maine be the data inwn: ae 7 torte Mean sentigieteyaaetinmcnysbnl etme 


sdolehbia x68 4 
(io) ot 

: 

EVEF = 
ref 
if ~~ 

Py hea 


. 
! 
: 
t- 
: 
: 
: 
: 
; 
} 
| 
in 0, ,8&- 
7. 


Di f her 


: 
oe 


se a al i le gis a 


7 | ; 
] 


(Baws pngen ser 
toy ok oat or ee v2 


a 4 i 0a Tae hy, ee | ee ee 
nt : i ral oh mie Vy a, maat a ny) re / 5 : 
j 1 ; ia} 
= eee Miki Yee : ; 


. 1 
at 
ac a 


= 


Per ; it be ‘ 
Cane 


AO . i : A 

ce ‘inte Ph Aa? 

Vid aia eee road 
eves ab 


ba Sas 


etsorial) aa  yexis: 36a rhe 
te 1 i Ay _ ; } A aH 
oF betasiee a va Ibiw wigs Sed 


sVitocxld sag Foo a sh ‘tend 


ne 


vA Ss 
ne date Getty ahs Bbotss . 


i Li The 


on -: 
rhe oo : $ 
Prahieir, . 


ais thee ‘ads od asa sat Bip 


= EL = 
VS ¥ 


d3hL"eesG | Si pat 
(seeTpad) [» (25oqped). 
| pe, ' 


eee 


> 
== 
uj a ie 


066.0. ts peso)” 1 


| eG 
bae | oF 20 .. Ms 
ra as : 
hie Poi ake al aati 
SOU bk eee 
a aa A Prorat 
BORO fy AGO, | 


fa, 


a = a we ae \ ve. e e ‘ 


seeks oA 
‘hot Gh nee 


beamwidth when the main beam is shifted up to 30 degrees 
from the zenith are shown in Figure 3.6. For @ close to 
zero the beamwidth of the Gaussian array is 46% greater than 
the beamwidth of the uniform array and for @ close to 30 
degrees it is 50% greater. The percentage increase in 
beamwidth is 12.8%. For the Gaussian array the percentage 


increase is 15.2%. 


3.4 Sidelobe Levels 

The purpose of tapering is to reduce the sidelobe 
levels of the array factor. Consider first the peak 
Sidelobe levels, as functions of 6, of the uniform array and 
the truncated Gaussian arrays shown in Figure 3.7. The 
curves represent the locus of values of the peaks of the 
Sidelobes. It is apparent that for both arrays the sidelobe 
levels decrease smoothly for increasing 6. In both cases 
they drop off rapidly from the level of their first sidelobe 
and then level out as @ approaches 90 degrees. The peak 
sidelobe levels for small values of § are about 25 dB lower 
in the case of the array with the Gaussian taper. But, for 
6 greater than about 2.5 degrees the difference is about 
20dB. It was observed during the study, although it is not 
shown here, that for large values of @ the location of the 
Sidelobes and the zeros in the array factors for both array 
systems occur very close together. The spacing between the 
zeros for both systems for large values of 9 is determined 


by the overall dimensions of the array. 
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Phased Location of the Main Beam (degrees) 


Figure "3.6. 
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Attempts were made to obtain similar diagrams for 
the three coarse approximations to the Gaussian excitation 
distribution but it was found that the peak sidelobe levels 
do not follow a smooth curve. The scale of 0 in Figure 3.7 
is insufficient to show these variations in the sidelobe 
levels. The array factors were then drawn in sections to an 
expanded scale of 6. Figures 3.8a through 3.8c show the 
peak sidelobe levels that occur in the first 10 degrees of 
6. The sections from 2 degrees to 10 degrees can be 
considered typical for each system and only the analysis of 


these sections will be presented here. 


One of the most striking features that can be seen 
on these diagrams is that there are points, in all three 
systems with course tapering, where the peak sideloke levels 
are greater than the peak sidelobe levels of a uniform 
array. What has happened is that the tapering has 
eliminated some of the sidelobes and reduced many of the 
others. Some sidelobes have remained almost unchanged. The 
reason that they appear to be greater than the sidelobes of 
a uniform array is because each array factor was normalized 
with respect to its own maximum. If they had all been 
normalized with respect to the maximum of the uniform array, 
the array factors for the arrays with two, four and eight 
stages of tapering would have to be reduced by 2.3 dB, 4.2 dB 
and 4.8 dB respectively. If this were done some of the peak 


sidelobe levels of the arrays with coarse tapering would 
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just touch the smocth curve through the peak sidelobe levels 


of the array with a uniform excitation distribution. 


The frequency with which the peak sidelobe levels 
of the array with coarse tapering are greater than the 
levels of the array with uniform excitation were measured. 
It was found that as the number of stages of tapering was 
increased the frequency at which large sidelobes occurred 
was reduced. Hence, increasing the taper reduces more 


sidelobes. 


The average sidelobe level for each of the array 
systems was calculated. The results, which are given in 
Table 3.2, show that as the number of stages of tapering is 
increased the average sidelcbe level is reduced. 

Table 3.2. 


Average Sidelobe Levels close 
to the Main Beam. 


(So cr p aeee es PRR ST PN BEM OLS AS 
{ Taper | Average Sidelcbe | 

{ Level (dB) | 
{ { { 
|-------- 4 --_-------—_- 
{ Uniform | -3185 | 
| | ( 
{| 2-stage | =66.4 | 
\ | i 
| 4-stage | *=38<2 | 
I { i 
{ 8-stage | -42.1 i 
\ { { 
{ Gaussian | -51.4 . 
| i 
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CHAPTER 4 


The array factors of several one-dimensional 
arrayS were examined in Chapter 3. These also could be 
considered as two-dimensional arrays with constant array 
factors in their transverse planes. If a number of 
identical one-dimensional arrays are placed side by side to 
form a two-dimensional array then the resulting array factor 
will be the product of two orthogonal one-dimensional array 
factors as indicated in equation 2-52. If the one- 
dimensional arrays are not identical then the array must be 
kroken up into sub-arrays where the one-dimensional arrays 
used in the sub-arrays are identical. The array factor for 


this type of array may be found by using equation 2-56. 


The physical layout of the elements of all the 
arrays to be discussed in this chapter will follow either 
the arrangment of an 8 elewent wide rectangular array shown 
in Figure 4.1a or of an 8 element wide array with four 
stages of physical tapering shown in Figure 4.1b. It should 
fe noted that it is possible, by a suitable choice of the 
excitation distribution, to design several different array 
systems combining both physical and resistive tapering such 


that the array factors in the principle plane are identical 


for each system. 
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The arrangment of the elements and the excitation 
distribution of the arrays are symmetric about both X and Y 
so that it is only necessary to examine the array factor in 
cne quadrant of K-space. In order to achieve a reasonable 
degree of resolution and to insure that no large sidelobes 
are missed, it is necessary to compute the array factor at 
approximately 106 points in K-space. A very considerable 
amount of computing is obviously required. The final output 
is arrived at in two stages. One quadrant of K-space, 
defined in the region 0<K,<2m7/A and OSK,<27/1, is divided up 
into 32 incriments along the K, and Ky axeS aS Shown in 
Figure 4.2. In turn this gives 1024 equal sections. The 
peak value of the array factor in each section, which also 
represents the peak sidelobe level in that section, is 
found. The section of K-space representing real space is 
now divided into 120 blocks along the 8 and @ axes. For 6 
less than 30 degrees there are 60 of the latter blocks where 
each block represents a change in @ cf 5 degrees and a 
change in @ of 10 degrees. For 9 greater than 30 degrees 


there are also 60 blocks but in this case each block 


represents a change in both 6 and ¢ of 10 degrees. 


4.2 Excitation 
The array factor of four two-dimensional arrays 
will be examined in this chapter. The physical layout and 


excitations of the arrays are shown in Figures 4.3a through 


4.3d. For convenience they will be named as follows 
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Block Representation of K-space. 


Figure 4.2 
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Figure 4.3a. Arrangement of the Elements and the 
Excitation Data for the Uniform Array. 
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Figure 4.3b. Arrangement of the Elements and the 
Excitation Data for the Physically 
Tapered Array. 
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Figure 4.3c. Arrangement of the Elements and the 
Excitation Data for the Combination 
Array. 
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Figure 4.3d. Arrangement of the Elements and the 
Excitation Data for the Gaussian 


Array. 
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Uniform Array: The shape is rectangular and the 
excitation is uniform. 


Physically Tapered Array: The shape is physically tapered 
and the excitation is uniform. 


Combination Array: The shape is physically tapered 
and the excitation is resistively 
tapered. 

Gaussian Array: The shape is rectangular and the 


excitation is resistively tapered 


The shape of the array factors in the principle 
plane of the two-dimensicnal uniform and Gaussian arrays are 
respectively identical to those of the one-dimensional 
uniform and Gaussian arrays, aS can be seen by setting P, to 
zero in equation 2-52. The level of the two-dimensional 
array factors is 9 dB greater than the level of the one- 
dimensional array factors. This increase is the result of 
the increase, in the two-dimensional case, of the number of 
elements to 8 across the width of the array. The shape of 
the array factor in the principle plane of the two- 
dimensional physically tapered array is the same as that of 
the one-dimensional array with four stages of tapering. The 
shape of the array factor in the principle plane of the two- 
dimensional combinaticn array is the same as that of the 
one-dimensional Gaussian array. In both cases the levels of 
the array factors of the physically tapered and the 


combination arrays are 9 dB higher than the levels of the 


corresponding one-dimensional arrays. 
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The shape of the array factor in the transverse 
plane depends on the amount of tapering that exists along 
the transverse axis. For a single line of dipoles the array 
factor in the transverse plane is constant. For a two- 
dimensional array made up of identical one-dimensional 
arrays placed side by side, such as the two-dimensional 
uniform and Gaussian arrays, the effective excitation along 
the transverse axis is constant as shown in Figure 4.4a. 
The array factor in the transverse plane for both of these 
arrays is shown in Figure 4.5a. The effective excitations 
shown in Figures 4.4a through 4.4c are found by summing C,, 


over all n at each value of nm. 


The effective excitation along the transverse axis 
for the two-dimensional physically tapered and combination 
arrays are Shown in Figures 4.4b and 4.4c respectively. The 
array factors in the transverse plane for these arrays are 
shown in Figures 4.5b and 4.5c respectively. They show a 
significant increase in the beamwidth and decrease in the 
sidelobe levels compared to those of the uniform or Gaussian 
arrays. The precise value of the decrease in the sidelobe 
level is not important. It wiblybe shown inthe next 
section that the sidelobes closest to the transverse plane 
for both the physically tapered and the combination arrays 
are not necessarily the greatest fcr arbitrary values of 06 
greater than 30 degrees, however, they are im’general Lower 


than the sidelobe levels of the uniform or Gaussian array in 


the same region. 
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r ox 
a. Distance from the Centre (A) 
LAO 
2), A 0 x 2X 
b. Distance from the Centre ()) 
LEO 
2 r 0 r 2X 
Cre Distance from the Centre (A) 


Figure 4.4. Effective Excitation in the Transverse 
Plane. 
a) Uniform and Gaussian Arrays. 
b) Physically Tapered Array. 
c) Combination Array. 
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4.4 Peak Sidelobe Levels 

The peak sidelobe levels for the four arrays 
described in Section 4.2 are given in Tables 4.1 through 
4.4. For the purpose of making comparisons between the 
different array systems a quadrant of the sky will be 
divided into five regions as indicated in the tables. 
Region 1 (9<5, 0<$#<90) represents the area close to main 
beam. Regions 2 and 3 represent the areas close to the 
principle plane and the transverse plane (80<$<90, 5<@<90 
and @<10, 5<6<90) respectively. Region 4 (5<0<30, 10<¢<80) 
represents the area off the principle planes for §$ less than 
30 degrees and Region 5 (30<0<90, 10<#<80) represents the 


area off the principle planes for @ greater than 30 degrees. 


Region 1 Main Bean. 


Region 1 contains the main beam. In the tables 
the main beam has been suppressed and the entries represent, 
for each array, the largest sidelobe found in the region. 

In the case of the uniform and the Gaussian arrays this 
Sidelobe is the first sidelobe in the principle plane, 
however, for the physically tapered or the combination 
arrays the location of the largest sidelobe is unknown. The 
relationship between the array ee the principle 
plane of the two-dimensicnal arrays and the array factors of 
the one-dimensional arrays was indicated in the previous 
section. The analysis of the main beams of the one- 


dimensional arrays with uniform and Gaussian excitation 
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Table 4.1. 


10 


Spherical Coordinate @ (degrees) 


20 30 40 50 60 70 80 


Maximum Sidelobe Level: i305 dB. 


Peak Sidelobe Levels in dB Below the Main 
Uniform Array. No ELrors. 
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Table 4.2. 
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Spherical Coordinate @ (degrees) 


20 30 40 50 60 70 80 


Maximum Sidelobe Level: 25.5 @B. 


Peak Sidelobe Levels in dB Below the Main 
Physically Tapered Array. No Errors. 
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Table 4.3. 
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Spherical Coordinate @ (degrees) 


20 30 40 50 60 70 80 


Maximum Sidelobe Level: 34.9 dB. 


Peak Sidelobe Levels in dB Below the Main 
Combinaticn Array. No Errors. 
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Table 4.4. 
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Spherical Coordinate @ (degrees) 


20 30 40 50 60 710 80 


Maximum Sidelobe Level: 40.8 dB. 


Peak Sidelobe Levels in dB Below the Main 
Gaussian Array. No Errors. 
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distributions and the array with four stages of tapering, 
considered in Section 3.3, is applicable to the arrays 
discussed in this chapter. For points off the principle 
planes but still close to the main beam equations 2-52 and 
2-56 must be used to find the array factors. By examining 
the F(P,) terms in these equations it can be seen that for 
small values of P, the dependence of the array factor on BS 
is weak. It can be seen from Figures 4.5a through 4.5c 
that, for 6 less than 5 degrees, the effects of P, on the 
level of the main beam of the four two-dimensional arrays 
will differ by less than 1 dB. The Gaussian array has the 
lowest maximum sidelobe level in this region. The maximun 
Sidelobe level of the combination array is 6 dB higher close 
to the transverse plane. Nc further analysis of this region 


will be made. 


Region 2, Principle Plane. 

In the region close to the principle plane the 
Gaussian array has the best sidelobe level performance. The 
peak sidelobe levels decrease from -52.9 dB for 9 close to 5 
degrees to -69.4 dB for @ close to 90 degrees. The 
performance of the combination array is not as good. There 
is little difference for small values of 0; however, for @ 
close to 90 degrees the peak sidelobe level is -56.3 dB. 
This represents an increase of 13.1 dB over the level of the 
Gaussian array. The performance of the uniform and the 
physically tapered arrays are simular to each other. The 


sidelobe levels go from -33.5 dB for small values of @ to 
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-50.3 dB for large values of 9. This represents a serious 
loss in the performance compared to the Gaussian or the 
combination arrays. For small values of @ the sidelobe 
levels of the unifcrm or the physically tapered arrays are 
about 20 dB higher than those of either the Gaussian or the 
combination arrays. For large values of @ the increase in 
Sidelobe level is about 19 dB compared to the Gaussian array 


and about 6 dB compared to the combination array. 


The physically tapered array has sidelobe levels 
up to 3.6 dB higher than those of the uniform array. As 
pointed out in Section 3.4 these sidelobe levels only appear 
higher than those cf the uniform array because of the way in 


which the respective array factors were normalized. 


Region 3 Transverse Plane. 

The peak sidelobe level characteristics near the 
transverse plane for the two-dimensicnal arrays fall into 
two catagories. For arrays with no physical tapering, such 
as the uniform cr the Gaussian array, the peak levels are 
almost identical. They also follow closely the array factor 
in the transverse plane, Figure 4.5a. This result may be 
anticipated since the array factors for these arrays may be 
expressed as the product of two orthcgonal array factors 
which both show a progressive decrease in the sidelobe 


levels for increasing 9. The peak levels decrease from -2.1 


dB for small values of 9 to -26.4 dB for large values of 6. 
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For the arrays with physical tapering, such as the 
physically tapered array or the combination array, the peak 
sidelobe levels, for ® greater than about 20 degrees, are 
greater than the values of the array factor plotted in the 
transverse plane, Figures 4.5b and 4.5c, over the 
corresponding ranges of @. This result indicates that for 
the physically tapered arrays the largest sidelobes in 
region 3 do not occur in the transverse plane. In general 
their locations are not easily predictable. Tables 4.2 and 
4.3 for these arrays indicate that there is a general 
decrease in the peak sidelobe levels from -1.3 dB for small 


values of 9 to -33.4 dB for large values of 6. 


For @ less than 20 degrees the peak sidelobe 
levels of the combination array are higher than those of the 
Gaussian or the uniform arrays. This result is due to the 
increase in beamwidth in the transverse plane that 
accompanies physical tapering. The maximum difference of 
4.0 GB occurs for @ between 10 and 15 degrees. For 9 
greater than 20 degrees the combination array produces the 
most rapid decrease in the peak sidelobe levels. The peak 
sidelobe levels of the combination array, for large values 


of 6, are about 7.0 dB lower than those of the Gaussian 


array. 


Region 4, Near Sidelobe Levels off the Principle Planes. 


In this region the Gaussian array has the lowest 


sidelobe levels. They vary from -45.1 dB to -72.3 dB. The 
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average level of the peaks is about -57.5 dB. For the 
combination array the peak sidelobe levels vary from -30.9 
dB to -53.5 dB with the average peak level at -42.5 dB. The 
Sidelobe levels of the combination array decrease steadily 
for increasing $9. The Gaussian array displays a similar 
decrease in the sidelobe levels but in addition the first 
zero of the array factor in the transverse plane introduces 
a large drop in the levels near 9 equal to 20 degrees and 9 


equal to 35 degrees. 


The peak sidelcbe levels of the physicaily tapered 
and the uniform arrays are higher than those of the Gaussian 
or the combination arrays. For $@ between 10 and 20 degrees 
there is little difference between the uniform or the 
physically tapered array and the combination array. For 9 
close to 80 degrees the peak sidelobe levels of the uniform 
or the physically tapered arrays are about 15.0 dB higher 


than those of the combination array. 


Region 5, Far Sidelobe Levels off the Principle Planes. 

In this region the Gaussian array has the lowest 
sidelobe levels. Over most of the region they are about 20 
dB lower than those of the uniform array and about 30 dB 
lower than those of either the physically tapered or the 


combination arrays. They vary from -62.7 dB to -85.8 dB. 


4.5 Discussion 
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Considering the main beam, the array factors in 


the principle planes and the peak sidelobe levels the 
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Gaussian array has the most satisfactory power pattern 
characteristics. The highest sidelobes that may be expected 
will be about 38 dB lower than the main beam. The 
performance of the combination array in the main beam area 
is almost identical to that of the Gaussian array. In the 
area close to the main beam the highest sidelobe levels that 
may be expected for the combination array are about 6 dB 
higher than those of the Gaussian array. Both the uniforn 
and the physically tapered arrays have sidelobes between 15 
and 25 dB higher than those of the Gaussian array. In the 
regions far from the main beam the Gaussian array has 
sidelobe levels 30 dB below the peak sidelobe levels of the 
combination or the physically tapered arrays and 20 dB lower 
than for the uniform array. It is in this region that 
random errors in the excitation of the elements are expected 
to have the greatest effects. This will be examined in the 


next chapter. 
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CHAPTER’ 5 


EFFECTS OF RANDOM ERRORS 


5-1 Introduction 

The analysis of the one and two-dimensional arrays 
performed in the earlier chapters assumed that there would 
be no deviation in the construction or operation of the 
arrays from ideal arrays of identical elements equally 
Spaced and having specified excitations. This situation 
will not occur in practice. The effects of the mutual 
impedances on the excitation and the gain of the elements 
will not be considered in this thesis. There are, however, 


several other sources of errors which can easily be 


accounted for in the calculation of the array factors. 


An important source of errors is related to the 
construction and tuning of the arrays. The elements of the 
arrays may not be absolutely identical. There may be errors 
in the location of the elements or in the lengths of the 
cables joining them to the receiver. Location or cable 
length errors of a fraction of a percent will introduce 
phase errors of 1 degree or more in the excitation of the 
element. The processing equipment may introduce errors in 
the signals from the elements before they are correlated. 
It may be possible to compensate for or to correct some of 
these errors during the tuning of the array. However, the 
accuracy of the tuning will depend on the specifications of 
the equipment used. 
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Another source of errors is related to the 
physical enviroment in which the array and the processing 
equipment are placed. In the Northern Hemisphere ambient 
temperatures may vary from a low of -50 °F in the Winter to 
a high of about 100 °F in the Summer. In general an array 
would also be subjected to coatings of dust, water, snow and 
ice, all of which may affect its performance. Obviously, 
without a high degree of stability these enviromental 
extremes can produce significant errors in the excitation of 


an array. 


If the errors arising from the construction, 
tuning or stability of the array systems are completely 
random they can be accounted for when computing the array 
factor by introducing equivalent random errors in the 
Magnitude and phase of the excitation of the elements. It 
will be assumed for the analysis of the arrays in this 
chapter that the combined effect of all the errors can be 
accounted for by introducing magnitude errors in the 
excitation of the elements which will not exceed 0.5 dB and 
phase errors which will not exceed 5 degrees. The choice of 
these limits is based on the fact that such tolerances on 
the array system could be achieved in practice with careful 
design and construction. If the performance of the arrays 
is degraded to an unacceptable level by errors of these 
magnitudes it will be necessary to examine carefully the 
source of the errors and to attempt to reduce them. For a 


given system it would be possible to find the expected 
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errors by analysing the construction and tuning technigues 
and the specifications of the equipment used in the system 


er for tuning. 


S.2 Generation of the Random Errors 
The term in equation 2-38 representing the 
excitation of the elements is Cons. 2o eecount for the 


errors C may be redefined as follows 


mn 
Con = Amn(1t+e,, ) Exp{jlyn} (5-1) 


where €,, and ~,, are random variables within the limits 


defined by equations 5-2a and 5-2b. 


-0.060 < enn 


1A 


0.060 | (5-2a) 


1A 


The calculation of the array factor will be performed on an 
IBM 360 so ce,, and j,, can most conveniently be generated 
by using the subroutine RANDU in the IBM Scientific 
Subroutine Package [16]. In order to see more clearly the 
effects of the errors, the magnitude and phase errors were 
applied separately. The analysis showed that magnitude 
errors of 0.5 dB degrades the performance of the system to 
almost the same extent as phase errors of 5 degrees. The 
effects of both of these errors on the average sidelobe 
levels and the peak sidelobe levels are so nearly the same 
that only the results of the analysis with the phase errors 


will be presented here. 
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The array factors were recomputed using four 
different sets of random variables. There were some large 
gaps between the peak sidelobe levels for any one set of 
random errors so that it was difficult to predict the 
maximum levels that may be expected in these regions. By 
using four sets of random errors it was found that there 
were sufficient data points available to predict the peak 
sidelobe levels in all regicns of K-space. It was also 
possible with the additional data to assure that there were 


no trends in the random errors. 


5.3 Effects of Errors on the One 
The five one-dimensional arrays discussed in 
Chapter 3 are considered here with random phase errors of 
5 degrees in the excitation of the elements. The gain, 
beamwidth, the location of the first minimums and the peak 
sidelobe levels for the arrays with no excitation errors 


were found in Chapter 3 and will now be compared to those 


obtained when excitation errors are taken into account. 


For all the systems the effect of the errors on 
the gain was less than 0.1% and the changes in beamwidth 
were not detectable. The only sections of the main beams to 
be affected by the errors were those close to the first zero 
where the levels were more than 40 dB lower than the level 


of the main beam. 


To examine the effect of errors on the peak 


sidelobe levels it is necessary to consider the area close 
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to the main beam separately. Figure 5.1a shows the first 

10 degrees of the array factors for the uniform and the 
Gaussian arrays. The solid lines represent the peak 
sidelobe levels found when errors were introduced into the 
excitation of the elements. The broken lines which 
represent the peak sidelobe levels for the arrays without 
errors are included for reference. It is apparent that the 
effect of the errcrs on both systems has been to increase 
the sidelobe levels. In the case of the uniform array the 
difference is small for 6 less than 4 degrees and increaes 
Slowly beyond this point. In the case of the Gaussian array 
the increase starts at about 2 degrees and is much more 
pronounced. However, for the iatter system the sidelobe 
levels are much lower to start with and it is to be expected 


that the effect cf the errors is greater. 


The peak sidelobe levels in the first 10 degrees 
of the array factors for the one-dimensional arrays with 2, 
4 and 8 stages of tapering and with errors of 5 degrees in 
the phase of the excitations are shown in Figures 5.1b 
through 5.1d. For reference, the corresponding peak 
sidelobe levels for the arrays with no excitation errors are 
shown in each figure and the peak sidelobe levels for the 
uniform and the Gaussian arrays having no excitation errors 
are shown in all the figures. The effects of the errors on 
these arrays is the same as the effect of the errors on the 
uniform and the Gaussian arrays. The sidelobe levels have 


increased and the increase is greatest at points furthest 
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from the main beam. In all cases the sidelobes most 


affected by the errors are those which were low initially. 


It can be seen that the errors have very little 
effect on the magnitude or location of the largest 
sidelobes. The magnitude and location of these sidelobes 
depends on the sizes of the steps used for tapering and the 


points at which the steps occur. 


Figure 5.2a shows, for large vaiues of 0, the peak 
Sidelobe levels that may be expected for the one-dimensional 
uniform and Gaussian arrays when errors in the excitation 
are introduced. The peak sidelobe levels of the same arrays 


with no excitation errors are included for reference. 


It can be seen that the increase in the sidelobe 
levels due to the excitation error which was observed in 
Figure 5.1a continues for large values of 6. For both 
arrays the peak sidelobe levels fall towards a limiting 
value. The limit is -41 dB for the Gaussian array and 1 dB 
higher for the uniform array. When 90 is 8 degrees in the 
case of the Gaussian array and 30 degrees in the case of the 
uniform array the sidelobe levels have fallen to within 1 dB 


of their limits. 


Figure 5.2b shows, for large values of 0, the 
region within which peak sidelobe levels can be expected to 
occur for the arrays with 2, 4 and 8 stages of tapering when 


errors in the excitation are considered. The peak sidelobe 
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levels expected for these systems are simular. For any one 
system some sidelobes should be at least as high as the 
lower solid line and none should be greater than the upper 
solid line. It can be seen that these limits follow closely 
the peak sidelobe levels that may be expected for the 
uniform array when excitation errors are considered. The 
peak sidelobe levels for the uniform and Gaussian arrays 


With no excitation errors are included for reference. 


The average sidelobe level for the five one- 
dimensional arrays, without excitation errors and with 


excitation errors, were computed as follows 


T/2 
AF(€)dé@ 
AF(avg) = §(min1) (5-3) 


T/2 - @(mint) 


where @(min1) is the first minimum of the array factor. 

The results, which were normalized with respect to the 
levels of the main beams, are given in Table 5.1. It can be 
seen that the effect of the errors on the average sidelobe 
level is to increase the level. The increase is small for 
the uniform array (1.4 dB) and the increase is large for the 
Gaussian array (11.1 dB). An important result which can be 
seen from the table is that the spread in the average 
sidelobe levels between the different array systems is 
reduced when the excitation errors are considered. The 


array most affected by the errors is the Gaussian array. 
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Table 5.1. 
Average Sidelobe Levels 
Over all Sidelobes. 


iam (0: 8 on Ceerte eo fo TV 
{| Taper { No Errors | Errors | 
l | (dB) [| (dB) { 
{ { | { 
-—__—___—_-+-__-—-___—-++ 4 
{ Uniform) | -45.8 | -44 4 | 
{ | | i 
{| 2-stage | -48.3 | =45.9 | 
{ | | | 
{ 4-stage | =5097 | -47.0 { 
l { i \ 
{ 8-stage | -56.1 | -48.4 | 
i l i { 
{| Gaussian | =6 142 { #50 23 { 
{ l { i 
[oak a We a ee Jd 


The effect of excitation errors on the performance 
of the two-dimensional arrays discussed in Chapter 4 will 
now be considered. AS in Chapter 4 only the peak sidelobe 
levels will be examined. The sky will also be divided into 
the same five principle regions. The results of this 
analysis for errors cf 5 degrees in the phase of the 


excitations are given in Tables 5.2 through 5.5. 


Region 1. Main Bean. 
It was indicated in the last section that 
excitation errors had very little effect on characteristics 
of the main beam of one-dimensional arrays. It is clear 
from that discussion and the discussion of the main bean 
area in Section 4.4 that for the two-dimensional array the 


effects of excitation errors on the main beams are small. 


The maximum sidelobe levels of the Gaussian, the combination 
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Table 5.2. 
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Spherical Coordinate @ (degrees) 
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Maximum Sidelobe Level: 42.9 4B. 


Peak Sidelobe Levels in dB Below the Main 
Uniform Array. 5° Phase Errors 
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Table 5.3. 
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Spherical Coordinate @ (degrees) 


20 30 40 50 60 70 80 90 


Maximum Sidelobe Level: 22.14 dB. 


Peak Sidelobe Levels in dB Below the Main Beam. 
Physically Tapered Array. 5° Phase ELrors. 
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Table 5.4. 
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Spherical Coordinate % (degrees) 


20 30 40 50 60 70 80 90 


Maximum Sidelobe Level: 2t.2 dB. 


Peak Sidelobe Levels in dB Below the Main Bean. 
Combination Array. 5° Phase Errors. 
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Spherical Coordinate ® (degrees) 


10 


Spherical Coordinate @ (degrees) 


20 30 40 50 60 70 80 


Maximum Sidelobe Level: 36.4 dB. 


10 


15 


20 


25 


30 


40 


50 


60 


70 


80 


90 
Tabie 5.5. 


Peak Sidelobe Levels in dB Below the Main 
Gaussian Array. 5° Phase Errors. 


90 


Beam. 
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and the physically tapered arrays have increased by abcut 3 


Principle Plane. 
In the region close to the principle plane the 
combination array has the lowest sidelobe levels, which vary 
from -39.9 dB to -43.9 dB. The sidelobe levels of the 
physically tapered array are about 4 dB higher than those of 
the combination array over most of this region, hovever, for 
@ less than 10 degrees they are 10.6 dB higher. The peak 
Sidelobe levels of the Gaussian array vary from -33.1 dB to 
-38.0 dB which indicates that the errors have introduced a 


serious degradation in the performance of this array. 


Region 3. Transverse Plane. 
The effect of the errors on the sidelobe levels 
close to the transverse plane for all the two-dimensional 


arrays is less than 0.1 dB. 


Region 4 Near Sidelobes off the Principle Planes. 

Over most of this region the performance of the 
Gaussian and combination arrays are similar. The largest 
sidelobe for the Gaussian array are at -34.0 dB and for the 
combination array are at -30.7 dB. ‘The peak sidelobe levels 
cf both the physically tapered array and the uniform array 
are larger than those of either the Gaussian array or the 


combination array. This is mest clearly evident close to 


the main beam. 
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Region 5. Far Sidelobes off the Principle Planes. 

In this region the peak sidelobe levels of both 
the Gaussian array and the uniform array vary from about 
—31 dB to about -51 dB while those of both the physically 
tapered array and the combination array vary from about 
—37 dB to about -57 dB. For 8 close to 30 degrees the 


Gaussian array has the lowest sidelobe levels. 


5.5 Discussion 

The analysis of the effects of excitation error on 
the one-dimensional arrays showed that they are small close 
to the main beam but that they increase for large values of 
8. For €@ close to 90 degrees the peak sidelobe level 
approaches a limiting value. It can be seen that the same 
trend is apparent for the two-dimensional arrays; however, 
Since the array factor in the transverse plane is not 
affected by the errors, the first zero of these array 
factors influence the sidelobe levels over much of K-space. 


This situation accounts for the lowest sidelobe levels in 


Tables 5.2 through 5.5. 


Considering the overall performance of the two- 
dimensional arrays the Gaussian array is best in the 
vicinity of the main beam. The combination array has 
sidelobes about 6 dB higher close to the main bean. 

However, close to the principle plane and for @ greater than 
5 degrees the perfcrmance of this array actually becomes 


better than the Gaussian array. This statement is also true 
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for the region off the principle planes and where @ is less 
than 30 degrees. Both the uniform array and the physically 
tapered array have large sidelobes close to the main bean 
(at least 15 dB higher than the Gaussian array). fThe 
uniform array has the highest sidelobes in all regions where 
@ is less than 30 degrees. The physically tapered array has 
the highest sidelobes in the regions off the principle 


planes and far from the main bean. 
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CHAPTER 6 


SUMMARY AND CONCLUSIONS 


In this thesis it was shown that the power pattern 
of an antenna array can be considered as the product of an 
accray factor and a reference or element pattern. The array 
factor accounts for the arrangment and excitation of the 
elements in the array. The reference pattern depends on the 
shape of the elements used in the array and the position of 


the reflecting screen. 


A technigue was developed for computing the array 
factor of one and two-dimensicnal arrays which is based on 
the use of the fast Fourier transform. Because the number 
of steps involved in the calculations using this type of 
algorithm is small it is faster and more accurate than the 


usual summing techniques. 


The technigue was applied to the analysis of one- 
dimensional arrays without excitation errors. The first 
array to be examined had a uniform excitation distribution. 
Arrays with tapering in 2, 4 and 8 stages were considered. 
These represent coarse approximations to a Gaussian 
distribution. Finally, an array with a Gaussian excitation 
was examined. As the excitation distribution was changed 
from uniform through progressively better approximations to 
a Gaussian, the following changes in the array factor were 


observed. 
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1. the gain decreased. 

2. the beamwidth increased. 

3. the average sidelobe level decreased. 

4. the peak sidelobe levels decreased. 

5. The arrays with 2, 4 and 8 stages of tapering had 
Sidelobe levels as much as 3 dB greater than those 
of the uniform array, however, the frequency with 
which these occurred decreased as the taper became 


more Gaussian. 


Extension of the analysis to a two-dimensional 
uniform array and three progressively better approximation 
to a two-dimensional Gaussian array showed very little 
change in the array factors in the principle plane. The 
sidelobe levels in the transverse plane of both the 
physically tapered array and the combination array decreased 
compared to those of both the uniform and Gaussian arrays. 
This decrease in the sidelobe levels is due to the fact that 
the introduction of physical tapering in two-dimensional 


arrays introduces a taper along the transverse axis. 


In the area close to the main beam the sidelobe 
performance of the Gaussian array is best. The peak 
sidelobe level of this array is at least 6 dB lower than 
that of the combinaticn array and at least 15 dB lower than 


that of both the uniform and physically tapered arrays. 


In the area far from the main beam the sidelobe 


level performance of the combination array is best. The 
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principle contribution to the power from this region comes 
from the area close to the transverse plane where the peak 
Sidelobe levels of the combination array are lower than 
those of the Gaussian array by 3 dB when 90 is 30 degrees and 


by 6 dB when 6 is 90 degrees. 


The effects of introducing magnitude and phase 
errors in the excitation of the elements were examined. In 
the case of the one-dimensional arrays it was found that 
magnitude errors of 0.5 dB and phase errors of 5 degrees 
produced comparable degradation in the performance of each 
system. The principle effect of the errors was to increase 
the average sidelobe level and the peak sidelobe levels. 

The effects of the errors on the main beam and the areas 
close to the main kteam were small. The array most affected 
by the errors was the Gaussian array which initialiy had the 
lowest sidelobes levels. The most Significient result to be 
cbtained from the analysis was the reduction in the spread 
in the average sidelote levels between the systems from 25.4 


dB to 6 dB by the effects of the errors. 


The introducticn of excitation errors in the two- 
dimensional arrays produced severe degradation in the 
performance of all the arrays. As was the case for the one- 
dimensional arrays the effect of the errors in the area 
close to the main beam was small. In this region the 
Gaussian array is still the best. When the regions far from 


the main beam were examined it was found that the most 
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Significient contribution to the power came from the area 
close to the transverse plane where the performance of the 
combination array is up to 6 dB better than that of the 
Gaussian array. The sidelobe level performance of both the 
Gaussian array and the combination array was found to be 
very much better than that of the uniform array or the 


physically tapered array. 


On the basis of these results it is concluded that 
the performance of an array combining both physical and 
resistive tapering is comparable to that of an array with 
only resistive tapering. The performance of an array with 
only physical tapering would not be satisfactory due to the 
large sidelobes close to the main beam which would be 


introduced by such a taper. 


If the arrays discussed in this thesis are to be 
used as part of a T-interferometer it would be necessary to 
consider the combined operation of both arms of the T. fhe 
sidelobes of the T cannot be obtained directly from the 
results presented in Chapter 5 as only the peak sidelobe 
levels were considered. If any further investigation of 
these types of arrays is to be conducted it would be useful 
to compute the average sidelobe level for the arrays in 


addition to the peak sidelobe levels. 
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APPENDIX A 


Referring to figure 2.1 the expression for |I-r'] 


is obtained as follows 


Xx = rSinéCos¢ 


Y= rSinesing 


Z = rcos®@ 
—> 
Behe = CANS to fy h) 
[c-r' (2 = (X' - rSin@Cosg@)2 + (¥' - rSinOSing)2 + r2Cosz¢@ 


= X'2 + yt2 - 2rSin@(X'Cos@ + Y'Sing) 


+ r2Sin2@(Cos29 + Sin?) + r2Cos20 
= r2 - 2rSinO(X'Cosd + Y'Sing) + r2 


= Del, — 25100 (X*Cos@: + Y Sind) (Ey 


i a 


In the far field 


(E"P << 1 - 2S5in0(X'Cosd + Y'Sin®g) 
i. 


jr-r' 2 = r2{1 - 2Sin@(X'Cos@ + Y'Sing) } 
r 


And 


" 


jf-F"| = c{1 - SinO(X'Cosg + Y'Sing)} 


iG 
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AFPENDIX B 


Equation 2-18 is evaluated as follows 


where 


and 


U(y'-Y!) 


co Cc C&C 


fof fo bToCuad (h! x8) 0 (ets 1!)6 (24) sin (ky 


-0CO —0O 


EXp{jK,X'+jK, ¥"} Jaxtaytaze 


L/2-(¥'-¥!| 


u(Y'-¥!'-L/2) - u(¥'-¥!+1/2) 


Integrating with respect to X' and Z* gives 


Aig (Rae Ky) 


where 
F 
Let Wt 
then awt 
F 

where 
G (L) 


i 


i] 


i Gin H< Pik eee 


if DG NT Sore (L7/2-|2"- VA) Exp iiky 184 aye 


fo (a (W'-1/2) -u (W'+L/2) } Sin {k (L/2-|W"{)} 


Exp {jKy (W'+Y t3aw' 


‘Exp {Jky ¥4}G (L) 


L72 
[Sin {k (L/2-(W"{) }Exp {jKyW'} aw" 
-L/2 


This is a standard integral evaluated in most elementry 
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texts on antenna theory [17]. Noting that 


K, = kSin6Cos¢ 


G(L) is given as 


G(L) = 2[Cos{(kL/2) SinOSin@} - Cos{kL/2} ] 
1 - Sin?6 Sinegd 


A (K 


mn K 


x! now becomes 


y) 


Aon (h Ky) = IoC yp ExptjK, 11+ IR, Y4} 6 (L) 


tt 


Cron EXP {GK X At IK YA} Ay (Ke Ky) 


where 


AL, (Ky 7 Ky) 21, (Cos {(kL/2) Sin@Sing} ~9Cos EkKL/ 2} 4 


=~" Sin Sin 
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APPENDIX C 


COMPUTER PROGRAMMS 


The following is a brief description of the 
computer programmes used to calculate the array factors of 
the arrays described in this thesis. A listing of the 
programmes 1S contained at the end of the Appendix. 
Subroutines which use standard techniques have not heen 


listed. 


MAIN PROGRANNE 


too 


This programme was used to perform the 
calculations involved in the analysis of the one-dimensional 


arrays. The input for the programme is 


1 iene & ccntrol parameter. 

2. ERRORM maximum magnitude error. 

3. ERROR Maximum phase error 

b. DL the points at which course tapering occurs. 


The first of the programme up to and including the 
section for applying the B(K) correction was used to 
generate the tables of transforms required for the analysis 


of the two-dimensicnal arrays. 


This programme was used to compute the array 
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factors and the peak sidelobe levels of the two-dimensional 
arrays. The input for the programme is 
1. Tables of transforms from the first main programme. 


2. Tables of COS functions. 
MAIN PROGRAMME 3 


This programme was used to find the peak sidelobe 


levels of the array factors as described in Section 4.1. 


AFACTS This package is an assembler subroutine used to 
compute the array factors and the peak sidelobe levels for 
the two-dimensional arrays. This portion was found to be the 
most time consuming part of all the caiculations. By writing 
the programme in assembler language the total number of 
machine steps required for the calculation is reduced to the 


minimum. 


CSL This subroutine is used to generate the tables for the 
SIN and COS of large angles which occur in the evaluation of 


equations 2-51 and 2-54. 


arrays. 


HARM This is the fast Fourier transform programme available 
in the IBM Scientific Subroutine Library. It is used in the 


evaluation of equations 2-51 and 2-54. 
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INTR1 This subroutine is used to perform the integrations 
needed for evaluating the average sidelobe levels described 


by equation 5-3. The trapezoidal rule is used. 


except that the integration is performed over the peaks of 
the sidelobe levels. It is less accurate but faster than 


INTR1. 


LOG10 This subroutine converts the array factor, which is 


normalized with respect to the peak level, from the 


magnitude form to the dB form. 


MAG1 This subroutine converts the array factor from the real 


and imaginary number form to the magnitude and phase form. 


MAX This subroutine is used to find the maximum value of the 


array factor. 


MINA1 This subroutine finds the location of the first 


Minimum of an array factor. 


LOCMAX This subroutine finds the location and magnitude of 


the all the sidelobes of an array factor. 


PITAB This subroutine generates a table of equally spaced 


numbers in the range -27 to 2T. 


RNUM This subroutine generates a table of random numbers 


with a uniform distribution. 
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RANM This subroutine is used to apply random errors to the 


phase of the excitation of the elements of 


described by equaticn 5-1b. 


an array as 


RANM This subroutine is used to apply random errors to the 


magnitude ofthe excitation of the elements 


described in section 5.2. 


of an array as 


phase of an array factor calculated using the subroutine 


HARM. 


SCALE This subroutine is used to normalize 


with respect to the maximum value found by 


SETA This subroutine is used to expand the 


so that the resolution of the array factor 


ZERO This subroutine is used to zero parts 


between programme steps. 


an array factor 


LOCMAX. 


Size of an array 


is increased. 


of storage 
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MAIN PROGRAMME 1 


DIMENSION A1(4096) ,A2 (4096), A3 (4096) ,A4 (4096), 
A5 (4096) , AREF (320) ,AR(320) , THETAT (1025), 
PIP (2048) ,COSPF (2048) ,SINPF (2048), 
M(3) ,INV (512) ,S (512) ,L (8), 
B1(400) ,B2 (400) ,B3 (400) ,B4 (400) ,B5 (400), 
L1 (400) ,L2 (400) , L3 (400) ,L4 (400) ,L5(400) 
EQUIVALENCE (B1(1),L1(1)),(B2(1) ,12(1)), (B3(1),L3(1)), 
(B4 (1) ,-L4(1)), (B5(1) -L5(1)) 
DATA N,NAI,NAR/2048,4096,320/, NL/8/,M/11,0,0/, 
PI, TAPER, D/3.1415927,-12.0, 0.5/ 


PROGRAMME TO EXAMINE THE EFFECT OF RANDOM ERRORS IN 
THE PHASE OR MAGNITUDE OF THE EXCITATION OF THE 
ARRAY ELEMENTS. 


READ IN DATA 


READ (5,101) L 
READ(5,102) ERRORM 
READ(5,102) ERRORP 
READ(5,101) IC 


GENERATE TABLES AND CONSTANTS. 


M = ERRORM/100.0 
P = ERRORP 
I 


180.0/PI 
= 1.0/1024.0 
CALL PITAB(PIP,PI,N) 
CALL CSL(PIP,COSPF,SINPF,N) 
DO 1 I=1,1025 
THETAT(I) = ARSIN (FLOAT (I-1) *R) 


GENERATE THE EXCITATION DATA FOR THE ARRAYS 
WITHOUT ERRORS IN THE MAGNITUDE OR THE PHASE. 


DO 99 IK=1,5 
CALL GENARY (AREF,NAR,L, NL, TAPER, Ik) 
Ic = 1 


ZERO THE ARRAYS. 


CALL ZERO (A1,NAT) 
CALL ZERO (A2,NAI) 
CALL ZERO (A3,NAI) 
CALL ZERO (A4,NAT) 
CALL ZERO (A5,NAI) 


EXPAND INTO A LARGER VECTOR TO INCREASE THE 
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C *ex 
Cex 
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RESOLUTION 


CALL SETA(A1,AREF, NAI, NAR) 
CALL SETA (A2,AREF,NAI,NAR) 
CALL SETA(A3,AREF, NAI,NAR) 
CALL SETA (A4,AREF, NAI, NAR) 
CALL SETA(A5,AREF, NAI,NAR) 


DETERMINE WHICH TYPE OF ERROR IS TO BE EXAMINED. 
MAGNITUDE IF Ic=1 
PHASE IF Ic=2 
NO ERRORS IF Ic=3 


ELTON TAL ZO B90 ELC 
GENERATE A RANDOM ERROR IN THE MAGNITUDE. 


vy) omy ig 

DK = EM 

CALL RNUM(AR,1.0,DK,1X1,NAR) 
CALL RANM(A2,AR,NAI, NAR) 
CALL RNUM(AR,1.0,DK,1X1, NAR) 
CALL RANM(A3,AR,NAI, NAR) 
CALL RNUM(AR,1.0,DK,1X1, NAR) 
CALL RANM(A4,AR,NAI, NAR) 
CALL RNUM(AR,1.0,DK, 1X1, NAR) 
CALL RANM(A5,AR,NAI, NAR) 
WRITE (6,103) ERRORM 

GO TO 30 


GENERATE A RANDOM ERROR IN THE PHASE. 


TRL WATT 

DK = EP*PI/180.0 

CALL RNUM(AR,0.0,DK,1X1,NAR) 
CALL RANE(A2,AR,NAI, NAR) 
CALL RNUM (AR,0.0,DK,1X1, NAR) 
CALL RANP(A3,AR,NAI,NAR) 
CALL RNUM(AR,0.0,DK,1X1,NAR) 
CALL RANP(A4,AR,NAI, NAR) 
CALL RNUM(AR,0.0,DK,1X1, NAR) 
CALL RANP(A5,AR,NAI, NAR) 
WRITE (6,104) ERRORP 


PRINT OUT THE EXCITATION CF THE ARRAY ELEMENTS. 


WRITE (6,107) 

DO 60 I=864, 1185 

N2 = 2*I 

N1 = N2-1 

WRITE(6,108) I,N1,N2, 
A1(N1),A1(N2) ,A2(N1) ,A2 (N2) ,A3(N1), 
A3(N2) ,A4(N1),A4(N2) ,A5 (N1) ,A5 (N2) 
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FIND THE ARRAY FACTOR BY FOURIER TRANSFORMING THE 
EXCITATIONS. 


CALL 
CALL 
CALL 
CALL 
CALL 


HARM (A1,M,INV,S,1,1FERR) 
HARM (A2,M,INV,S,1,1FERR) 
HARM (A3,M,INV,S,1,1FERBR) 
HARM (A4,M,INV,S,1,IFERR) 
HARM (A5,M,INV,S,1,1FERR) 


APPLY THE B(K) CORRECTION. 


CALL 
CALL 
CALL 
CALL 
CALL 


FIND 


CALL 
CALL 
CALL 
CALL 
CALL 


FIND 


CALL 
CALL 
CALL 
CALL 
CALL 


ROTATE (A1,COSPF,SINPF,-1,NAI,N) 
ROTATE (A2,COSPF,SINPF,-1,NAI,N) 
ROTATE (A3,COSPF,SINPF,-1,NAI,N) 
ROTATE (A4,COSPF, SINPF,-1,NAI,N) 
ROTATE (A5,COSPF,SINPF,-1,NAI,N) 


THE MAGNITUDE AND PHASE OF THE ARRAY FACTOR. 


MAG1(A1,NAZ) 
MAG1(A2,NAT) 
MAG1 (A3,NAI) 
MAG1(A4, NAT) 
MAG1(A5,NAI) 


THE MAXIMUM VALUE CF THE ARRAY FACTOR. 


MAX (A1,NAI,A1MAX,LA1MAX, 2) 
MAX(A2,NAI,A2MAX,LA2MAX,2) 
MAX (A3, NAL, A3MAX, LA3MAX, 2) 
MAX (A4,NAI,A4MAX, LASMAX,2) 
MAX (A5,NAI,A5MAX, LASMAX, 2) 


NORMALISE THE ARRAY FACTOE. 


CALL 
CALL 
CALL 
CALL 
CALL 


FIND 


CALL 
CALL 
CALL 
CALL 
CALL 


FIND 


CALL 
CALL 


SCALE (A1,NAI,A1MAX, 2) 
SCALE(A2,NAI,A2MAX,2) 
SCALE (A3,NAI, A3MAX, 2) 
SCALE (A4,NAI,A4MAX,2) 
SCALE (A5,NAI, A5MAX, 2) 


THE FIRST MINIMUMS OF THE ARRAY FACTOR. 


MINA1 (A1,LA1MAX,LA1MIN, 2,N) 
MINA1(A2,LA2MAX,LA2MIN,2,N) 
MINA1 (A3, LA3MAX,LA3MIN, 2,N) 
MINA1(A4,LA4MAX, LAYMIN,2,N) 
MINA1 (A5, LASMAX, LASMIN, 2, N) 


THE PEAK SIDELOBE LEVELS. 


LOCMAX (A1,B1, 81,400, NEBL1,LA1MIN,N, 2) 
LOCMAX (A2,B2,B2,400, NEBL2,LA2MIN,N, 2) 
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CALL LOCMAX (A3,B3,B3,400, NEBL3,LA3MIN, N, 2) 
CALL LOCMAX(A4,B4,B4,400, NEBL4,LAUMIN,N,2) 
CALL LOCMAX(A5,B5,B5,400, NEBL5,LA5MIN, N, 2) 
XK OK 
*** FIND THE AVERAGE SIDELOBE LEVEL AS A FUNCTION OF 
eek THETA, 
aK 3K aK 
CALL INTTR3(A1,LA1MIN, THETAT,1025,SAT1) 
CALL INTTR3(A2,LA2MIN, THETAT, 1025,SAT2) 
CALL INTTR3(A3,LA3MIN, THETAT,1025,SAT3) 
CALL INTTR3 (A4,LA4MIN, THETAT, 1025,SAT4) 
CALL INTTR3(A5,LA5MIN, THETAT,1025,SAT5) 
SAT1 = SAT1/ (THETAT (1025) -THETAT (LA 1MIN/2+1) ) 
SAT2 = SAT2/ (THETAT (1025) -THETAT (LA2MIN/2+ 1) ) 
SAT3 = SAT3/(THETAT (1025) -THETAT (LA3MIN/2+1) ) 
SAT4 = SAT4/ (THETAT (1025) -THETAT (LA4MIN/2+1) ) 
SAT5 = SAT5/(THETAT (1025) -THETAT (LASMIN/2+1) ) 


%K KK 
*** FIND THE AVERAGE OF THE PEAK SIDELOBE LEVELS AS A 
*** FUNCTION OF THETE. 
KK IK 
CALL INTTR4 (B1,B1,THETAT, SBT1, NEBL1) 
CALL INTTR4(B2,B2,THETAT, SBT2,NEBL2) 
CALL INTTR4 (E£3,B3,THETAT, SBT3, NEBL3) 
CALL INTTR4(B4,B4, THETAT, SBT4,NEBL4) 
CALL INTTR4 (B5,B5,THETAT, SBT5, NEBLS) 
SBT1 = SBT1/ (THETAT (L1 (NEBL1- 1)) -THETAT (L1(1))) 
SBT2 = SBT2/(THETAT (L2 (NEBL2-1) ) -THETAT (L2 (1) )) 
SBT3 = SBT3/ (THETAT (L3 (NEBL3-1) ) -THETAT (L3(1))) 
1))) 
1))) 


SBT4 = SBI4/(THETAT(L4 (NEBL4-1)) -THETAT (L4 ( 
SBT5 = SBT5S/ (THETAT (L5 (NEBL5-1)) -THETAT (L5 ( 
aK 3K ak 
*** FOR EACH ARRAY PRINT OUT: 
a 2k THE LOCATION OF THE MAXIMUM LEVEL, 
kK THE VALUE OF THE MAXIMUM LEVEL, 
4K THE LOCATION OF THE FIRST MINIMUM, 
tok FIND THE AVERAGE OF THE PEAK SIDELOBE LEVELS AS A 
4K FUNCTION OF THETA. 
4K THE AVERAGE SIDELOBE LEVEL AS A FUNCTION OF THETA. 
aK 3K 3 
WRITE (6,110) 
WRITE(6,111) LAIMAX,A1MAX,LA1MIN,SAT1,SBT1 
WRITE (6,111) LA2MAX,A2MAX,LA2MIN, SAT2,SBT2 
WRITE(6,111) LA3MAX,A3MAX,LA3MIN,SAT3,SBT3 
WRITE (6,111) LA4SMAX, A4SMAX, LA4MIN, SAT4,SBTY 
WRITE(6,111) LASMAX,ASMAX,LA5MIN,SAT5,SBT5 
aK 3 3k 
*** WRITE OUT THE PEAK VALUES OF THE SIDELOBES. 
aK 3K 
WRITE(6, 120) 
DO 85 I=1,400,2 
J = I+1 
85 WRITE(6,121) B1(I) ,B1(d) ,B2 (I) ,B2(J) ,B3(I) ,B3(J), 
1 B4 (I) ,B4(J) ,B5 (I) ,B5 (J) 
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C *xKx& 
C *** CONVERT THE OUTPUT TO DB. 
C *x*x 
CALL LOG10(A1,NAI) 
CALL LOG10 (A2,NAI) 
CALL LOG10(A3,NAI) 
CALL LOG10 (AU, NAT) 
CALL LOG10(A5,NAI) 
C **x 
C *** FIND THE PEAK SIDELOBE LEVELS IN DB. 
C **x 
CALL LOCMAX(A1,B1,B1,400, NEBL1,LA1MIN, N,2) 
CALL LOCMAX(A2,B2,B2,400, NEBL2,LA2MIN,N,2) 
CALL LOCMAX (A3,B3,B3,400, NEBL3,LA3MIN, N, 2) 
CALL LOCMAX(A4,B4,B4,400,NEBL4,LA4MIN,N,2) 
CALL LOCMAX(A5,B5,B5,400, NEBL5,LA5MIN, N, 2) 
C *** 
C *** FOR EACH ARRAY PRINT OUT: 
C **x THE LOCATION OF THE MAXIMUM LEVEL, 
C #ex THE VALUE OF THE MAXIMUM LEVEL, 
C ex THE LOCATION OF THE FIRST MINIMUM, 
C *** 
WRITE(6, 110) 
WRITE (6,111) LAIMAX,A1MAX, LAIMIN 
WRITE(6,111) LA2MAX,A2MAX,LA2MIN 
WRITE (6,111) LA3MAX,A3MAX, LAZ3MIN 
WRITE(6,111) LA4YMAX,A4MAX, LA4MIN 
WRITE (6,111) LA5SMAX,AS5SMAX, LASMIN 
C **x* 
C *** WRITE OUT THE ARRAY FACTOR IN DB. 
C *x*x 
WRITE (6,112) 
DO 90 1=1025,1250 
{T = FLOAT (I-1025) 
J = IABS(I-1025) 
V = FLOAT (J) 
re 
IF (I.LE.1024) J = 141024 
AK = PI2*T*R 
THETAJ = ARSIN(V¥*R) *Q 
N2 = 2*J 
N1 = N2-1 
90 WRITE(6,113) I,J,AK, THETAJ, 
1 A1(N1),A1(N2) ,A2(N1) ,A2 (N2) ,A3(N1), 
2 A3(N2) ,A4(N1) ,A4(N2) ,A5 (N1),A5 (N2) 
C **x 
C *** WRITE OUT THE PEAK VALUES OF THE SIDELOBES IN DB. 
C *ex 
WRITE (6,122) 
DO 95 I=1,400,2 
J = I+1 


95 WRITE(6,121) B1(I),B1(dJ) ,B2(I) ,B2(J) ,B3(I) ,B3(J), 
1 B5 (I) -B5 (J) 
IF (IC.£0. 2). Go To 99 
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GO TO 5 
99 CONTINUE 
WRITE(6, 114) 
STOP 
C *** 


C *** FORMAT STATEMENTS. 
C **x 
101 FORMAT (8110) 
102 FORMAT (F10.6) 
103 FORMAT ("1EXAMINING THE EFFECT OF RANDOM MAGNITUDE!, 
1 ' ERRORS OF',F6.1, "PERCENT. ') 
104 FORMAT (' JEXAMINING THE EFFECT OF RANDOM PHASE ERRORS', 
' OF',F5,1,' DEGREES. *) 
cee FORMAT (* JEXCITATION DATA FOR THE ARRAYS. 


1 //,T12,*LOCATION',T28,"*TYPE 1 ARRAY',T50, 

2 'TYPE 2 ARRAY',T72,'TYPE 3 ARRAY',T94,'"TYPE 4', 
3 ' ARRAY',T116, 

4 'TYPE 5 ARRAY',//,' ELEMENT REAL IMAG',T27, 

5 ‘REAL IMAG',T49, "REAL IMAG! ,T71, 

6 ‘REAL IMAG', 

7 T93,"REAL IMAG',1T115, "REAL IMAG',//) 


108 FORMAT (16,19,15,5(F12.6,F10.6)) 
110 FORMAT (*1ARRAY FACTOR DATA.',////' LAMAX',T12,"AMAX', 
1 T26,*LAMIN',T41,*A(MEAN)',T61,"AP (MEAN) ', 
2 T81,*A(MEAN) ',1T101,"AP(MEAN)',//) . 
111 FORMAT (17, F14.7,110,4F20.10,/) 
112 FORMAT(*1ARRAY FACTORS FOR THE ARRAYS.', 
J//ei32,*TYPE 1 ARRAY',1T53,"TYPE 2 ARRAY',T74, 
'TYPE 3 ARRAY',T95,'TYPE 4 ARRAY',T116,'TYPE ¢, 
'S ARRAY',//,T5,%1',T10,'d",T14,"K',T19, 
‘THETA, 
730,*MAGNITUDE PHASE',1T51,'*MAGNITUDE PHASE', 
T72,'*MAGNITUDE PHASE',T93,'MAGNITUDE PHASE', 
T114,*MAGNITUDE PHASE',//) 
113 FORMAT (215,2F7.3,5(F12.3,F9.3)) 
114 FORMAT ("1") 
115 FORMAT (215, 2F7.3,5(F13.7,F8.3)) 
120 FORMAT ('1PEAK SIDELOBE LEVELS!,//) 
121 FORMAT (5 (110,F12.7)) 
122 FORMAT(*1PEAK SIDELOBE LEVELS IN DB.',//) 
123 FORMAT (5 (110,F13.3)) 
END 
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MAIN PROGRAMME 2 


DIMENSION A(9216) ,C40S6) ,AF (2048) ,AMAX (64) 


PROGRAMME TO COMPUTE THE ARRAY FACTOR OF TWO- 


DIMENSIONAL ARRAYS. 


A TRANSFORM TABLES. 
Cc COS TABLES. 
AF ARRAY FACTOR. 


AMAX PEAK LEVELS OF AF. 


EXTERNAL GETFD 

INTEGER*2 LEN 

INTEGER*4 FDUB,ADROF, FSIZE, BLOCK 
REAL*8 FNAME 

DATA NEB, NB/32,32/,LEN/8192/ 


READ IN THE TABLES OF DATA. 


READ(4) A 
READ(4) C 
READ(5,101) FNAME,FSIZE 


CREATE THE OUTPUT FILES. 
CALL GFILES (FNAME, FSIZE,FDUB,1,&501,6502) 
COMPUTE THE ARRAY FACTOR. 


READ(5,102) KSTART 

READ(5,102) NBLKS 

CM = C(1)*A (2) +C (2) *A (4) +C (3) ¥A (6) +C (4) *A (8) 
CP = C(1) *A(3) +C (2) ¥A (5) +C (3) *A (7) #C (4) *A (9) 
CNORM = 1.0/SQRT (CM*CM+CP*CP) 

K = KSTART 

CALL ZERO (AMAX ,64) 

BLOCK = (K-1)/32+1 

WRITE (10,112) 

DO 50 I=1,32 

KC = 4*K=3 

CALL AFACTS(A,AF,AMAX,C (KC) , NEB, NB,CNORM) 
K = Kt 

WRITE(9,111) AMAX 

DO 60 I=1,64,2 

IB = (I+1)/72 

I1 = I+1 

WRITE(10,113) BLOCK,IB,AMAX (I) ,AMAX(I1) 
NBLKS = NBLKS-1 

IF (NBLKS.GT.0) GO TO 1 

STOP 
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WRITE (6,601) 
STOP 64 
WRITE (6,602) 
STOP 68 


FORMAT STATEMENTS 


FORMAT (A8,18) 

FORMAT (10110) 

FORMAT ('1"',75,"B*,T10,'K',T27,"WRITE POINTER',/,"0',/) 
FORMAT (215,120) 

FORMAT (8X, 8Z8) } 

FORMAT('1 BIB ARMAX APMAX!,/,'0',/) 
FORMAT (215, 2F12.7) 

FORMAT(" RCALL FAILURE:"') 

FORMAT(' CREATE FAILURE:') 

END 
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MAIN PROGRAMME 3 


DIMENSION AB (2048) ,A1(512) ,A2(512) ,B1(512) ,B2(512) 


READ(5,101) AB 


IAB = 1 
I1= 1 
IL = 32 


pO 30 J=1,32 
DO 20 I=1I1,IL,2 


A1(IAB) = AB(TI) 
B1(IAB) = AB (I+1) 
A2(IAB) = AB (I+32) 
B2 (IAB) = AB(I+33) 
IAB = IABt1 

I1 = 11+64 

IL = IL+64 


WRITE(7,103) A1 
WRITE(7,103) A2 
WRITE (7,103) B1 
WRITE(7,103) B2 
CALL LOG10 (AB, 2048) 


IAB = 1 
I1 = 1 
IL = 32 


DO 50 J=1,32 
DO 40 I=11,IL,2 


A1(IAB) = -AB(I) 
B1(IAB) = -AB(I+1) 
A2 (IAB) = -AB(I+32) 
B2(IAB) = -AB(I+33) 
IAB = IAB+1 

I1 = 11+64 

IL = IL+64 


WRITE(7,102) A1 

WRITE(7,102) A2 

WRITE(7,102) B1 

WRITE(7,102) B2 

STOP 

FORMAT (8X,8Z8) 

FORMAT('1",/,© ©,16F5.1,(/,'0"s/e' 
FORMAT ('1*,16F8.5, (//,'0', 16F8.5) ) 
END 


1, 16F5.1)) 
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SUBROUTINE AFACTS 


CALL AFACTS (A, AF,AMAX,C,NEB, NE,CNORM) 


STM 
BALR 
USING 
ST 
LA 

LA 


14,12, 12(13) 


12,0 

*,12 

+3 ,SGR13 

13, SAVEAREA 

2,36 F*36! 
5,16(1) ADROF NEB 
5,0 (5) NEE 
5,=F*'1 NEB-1 

4,2 36* (NEB- 1) 
3,5 36% (NEB-1) 
4,8 F'gt 

5,20 (1) ADROF NB 
£,0(5) NB 

5,=F'1! NB-1 

5,3 8* (NB-1) 
6,0 (1) ADROF A 
3,6 ADROF A + 36*(NEB-1) 
7,8(1) ADROF AMAX 
6.7 ADROF AMAX + 8*(NB-1) 
8,12(1) ADROF C 
9,4 (1) ADROF AF 
10, 24(1) ADROF CNORM 
0,0 (6) AM 

0,0(8) AM*COS8 
0,0 AFRTM 
0,0(9) AFRTM 

0,0 (7) 

«+8 

0,0 (7) 

0,4 (6) A8SRE 

2,8 (6) ASIM 

0,0 (8) A8SRE*COS8 
2,0 (8) A8IM*COS8 
4,12(6) AGRE 

6,16 (6) A6IM 

4,4 (8) A6RE*COS6 
6,4 (8) A6IM*COS6 
0,4 

2,6 

4,20(6) AGRE 

6,24 (6) A4IM 

4,8 (6) AWRE*COS4 
6,8 (6) A4IM*COS4 
0,4 
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ESQRT 
ARG 


1,PARLSORT 
15, ESQRT 
14,15 
0,0(10) 
0,4 (9) 
0,4(7) 
*+8 
2,4(7) 
6,2,LOO0P 
3,1152 
7,8(7) 
7,4,LOOP 
13,SGR13 


14840940013) 


14 

F 

18F 

x' sor 
AL3 (ARG) 
V (SQRT) 
F 


A2RE 
A2IM 
A2RE*COS2 
A2IM*COS2 


LOAD PARLIST FOR SQRT 
LOAD ENTRY PT FOR SQRT 
BRANCH TO SQRT 
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SUBROUTINE GENARY(A,NA,L,NL, TAPER, TYPE) 
INTEGER*4 TYPE 

DIMENSION A(NA) ,L(NL) 

GO TO (10,20,30,40,50), TYPE 


INCJ = 8 
GO TO 60 
INCJ = 4 
GO TO 60 
INCJ = 2 
GO TO 60 
INCJ = 1 
GO TO 60 


Cc 


C = (~ALOG(10.0** (TAPER/10.0))) 
= C/(157.0*157.0) 


DO 51 I=2,160 


A(160+1) 
A (161-1) 
A(160) = 
A(161) = 
RETURN 
NI1 = 1 


EXP (-C* (I-1) * (I-1) ) 
A (160+TI) 


DO 280- dS 1 pO ZINC 

NIL = L(J+INCJ-1) 

C = 1.0-0.125* FLOAT (J-1) 
DO 70 I=NI1,NIL 


A (160+) 
A (161-1) 


C 
Cc 


NI1 = NIL+1 


RETURN 
END 
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SUBROUTINE CSL(A,COSL,SINL,N) 
DIMENSION A(N) ,COSL(N) ,SINL(N) 
Cc = 1.0 

DO 10 I=1,N 

COSL(I) = C¥#COS(A(I)) 

SINL(I) = -C¥SIN(A(I)) 

Gxstc¢ 

RETURN 

END 


SUBROUTINE ROTATE (A,COSPF,SINPF,IC,N,NCS) 
DIMENSION A(N) ,COSPF (NCS) ,SINPF(NCS) 
IF(IC.EQ.-1) GO TO 20 

DO 10 I=1,NCS 

X = A(2*I-1) *COSPF (I) A (2*1) *SINPF (I) 

Y = A(2*I) *COSPF(I) + A(2¥*I-1) *SINPF (1) 
A (2*I-1) = X 


A(2*I) = Y 
RETURN 
DO 30 I=1,NCS 


X = A(2*I-1)*COSPF(I) + A(2*1I) *SINPF (I) 
Y = A(2*1I) *COSPF (I) 
A(2*I-1) = X 

B(QeT yi ae ny 

RETURN 

END 


A (2*I-1) *SINPF (1) 
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SUBROUTINE RNUM(A,C,D,1X,N) 
DIMENSION A(N) 

DO 10 I=1,N 

CALL RANDU(IX,1IY,YFL) 

A(I) = (YFI-0.5) *2.0*D+C 
IX = IY 

RETURN 

END 


SUBROUTINE RANM(A,R, NA, NR) 
DIMENSION A(NA) ,R (NR) 

N1 = NA/2-NR+1 

DO 10 I=1,NR 

A(N1) = A(N1) *R (1) 


N1 = N1+2 
RETURN 
END 


SUBROUTINE RANP (A,R,NA, NR) 
DIMENSION A(NA) ,R(NR) 
N1 = NA/2-NR+1 


N2 = N1+1 

DO 10 I=1,NR 

C = COS(R(I)) 

S = SIN(R(I)) 

X = A(N1) *C-A(N2) *S 


Y = A(N1) *S+A(N2) *C 
RQ) =: xX 

A(N2) = Y 

N1 = N1+2 

N2 = N2+2 

RETURN 

END 
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